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Abstract. We consider relativistic many-particle operators which - accord- 
ing to Brown and Ravenhall — describe the electronic states of heavy atoms. 
Their ground state energy is investigated in the limit of large nuclear charge 
and velocity of light. We show that the leading quasi-classical behavior given 
by the Thomas-Fermi theory is raised by a subleading correction, the Scott 
correction. Our result is valid for the maximal range of coupling constants, 
including the critical one. As a technical tool, a Sobolev-Gagliardo-Nirenberg- 
type inequality is established for the critical atomic Brown-Ravenhall operator. 
Moreover, we prove sharp upper and lower bound on the eigenvalues of the 
hydrogenic Brown-Ravenhall operator up to and including the critical coupling 
constant. 
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1. Introduction and main result 



The description of atoms and molecules, in particular of their energies, has been 
a primer for the development of quantum mechanics. However, it became soon 
clear that atoms with more than one electron are not accessible to explicit solutions. 
This motivated the development of approximate models for large Coulomb systems. 
One of the most simple and - simultaneously - the most fundamental models was 
introduced by Thomas [S^ , Fermi [531 US] i and Lenz [3S] who proposed the energy 
functional which we will also use here. It predicts that the ground state energy 
of atoms would decrease with the atomic number Z to leading order as Z'' 1'^ . In 



order to get a refined description, Scott |49] conjectured that the electrons close 

to the nucleus should raise the energy by Considerably later Schwinger 

[47] argued also for Scott's prediction; Schwinger [48] and Englert and Schwinger 
[TPl [TTl [T^ even refined these considerations by adding more lower order terms 
[48] (see also Englert [9]). The challenge to address the question whether the 
predicted formulae would yield asymptotically correct results when compared with 
the iV-particle Schrodinger theory was for a long time unsuccessful. It was Lieb 
and Simon who proved in their seminal paper |39| that the prediction of Thomas, 
Fermi, and Lenz is indeed asymptotically correct. However, establishing the Scott 
correction resisted the mathematical efforts and became Problem lOB of Simon's 
15 Problems in Mathematical Physics [57]. Eventually, the Scott correction was 
established mathematically by Hughes [33l [34] (lower bound), and Siedentop and 
Weikard [50l [SH [52l [53l [54] (lower and upper bound) . In fact even the existence of 
the 2-5/3 -correction conjectured by Schwinger was proved by Fefferman and Seco 
[20l[ai[22l[l5l[23l[18l[l6l[17l[19]. Later these results were extended in various 
ways, e.g., to ions and molecules. 

Despite of the mathematical success in establishing the large Z asymptotics of 
the Schrodinger theory, these considerations remain questionable from a physical 
point of view, since large atoms force electrons into orbits that are close to the nu- 
cleus where the electrons move with high speed which should require a relativistic 
treatment. The atom is shrinking with increasing Z: already in non-relativistic 
quantum mechanincs the bulk of the electrons has a distance Z~^l^ from the nu- 
cleus; the electrons contributing to the Scott correction even live on the scale Z~^ . 
Schwinger [48] has estimated these effects concluding that a correction to the Scott 
correction occurs whereas the leading term should be unaffected by the change of 
model. S0rensen [46] was the first who proved that the latter is indeed the case 
for a simplified ad hoc naive relativistic model, the Chandrasekhar multi-particle 
operator, in the limit of large Z and large velocity of light c. In a previous paper 
[77] we established the value of the Scott correction which is again of order Z^, a 
result which was independently announced by Solovej, S0rensen, and Spitzer [59] 
(see also S0rensen [45] for the non-interacting case). Nevertheless, a question from 
the physical point of view remains: Although the Chandrasekhar model is believed 
to represent some qualitative features of relativistic systems, there is no reason to 
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assume that it should give quantitative correct resuhs. Therefore, to obtain not 
only qualitatively correct results it is interesting, in fact mandatory, to consider a 
Hamiltonian which ~ as the one by Brown and Ravenhall [4' - is derived from QED 
such that it yields the leading relativistic effects in a quantitative correct manner. 
(See also Sucher [62l[63l[64].) The first step in this direction was taken by Cassanas 
and Siedentop [5] who showed that, similarly to the Chandrasekhar case, the lead- 
ing energy is not affected. To show in which way the Scott correction is changed 
for this model is our concern in this paper. 

1.1. Relativistic energy form. According to Brown and Ravenhall [4] the energy 
of an atom with TV electrons in a state tp e 0^ is given by 




N 

^(cq;^-P^ + c2/3^-c2-Z|x^|-i)+ |x^-x^|"1 

1^=1 l<fj.<iy<N 



This involves the free Dirac operator reduced by the rest mass, acting in L'^(M.^, C*) 
with the four Dirac matrices in standard representation, 

"^(^ o)' -1 

where a are the three Pauli matrices in standard representation, i.e., 
1\ fO -i\ (\ 



'"1 = U oj' '"2= l,i oj' '"^ vo -1 

We use atomic units in which m — e? — h — \. The parameter Z is the atomic 
number and c the velocity of light. 

The Hilbert space of an electron is chosen as the positive spectral subspace of 
the Dirac operator, 

io"" := X[c^co)(ca • p + c2/3) (L2(M^ C^)) , 

and, correspondingly, the Hilbert space of N electrons f)^ is the antisymmetric 
tensor product of the one-particle space, i.e., f)^ :— A^Li-^l^- Fhially, the form 
domain of (P) is := n S(M3^,C^") with e the Schwartz space of rapidly 
decreasing functions. As is shown in [Tl], the Brown- Ravenhall form is closable 
and bounded from below if and only if 

2 k:=-<k^:=— -. 

^ ' c ~ 2/7r + 7r/2 

(See also Tlx [S71 [55] who improved the bound given in ^ to an explicite positive 
bound.) For the physical value, about 1/137, of the Sommerfeld fine structure, 
which equals 1/c in atomic units used here, the critical atomic number Z exceeds 
124 slightly. This includes all known elements. 

In the following we will assume that the atom described by H]) is neutral, i.e., 
Z — N, an assumption that we make mainly for the sake of brevity and clarity 
of presentation, since the Scott correction is independent of the ionization degree 
N/Z > const > 0. Similarly, it might seem that our treatment is restricted to 
spherically symmetric systems (atoms). However, on the energy scale considered 
here, molecular Hamiltonians essentially separate - in nature the distances between 
nuclei with charges ZZi, ZZk remain on a scale much larger than Z~^^^ ~ into 
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spherically symmetric one-center problems (atoms). Therefore, the molecular case 
follows from the atomic case by additional localization. However, for the sake of 
brevity and clarity, we will spare the reader the corresponding tedious technicalities, 
restrict to the atomic case, and freely use the resulting symmetry. 

Thus, according to Friedrichs, the one-particle form £^ defines for k < a 
distinguished self-adjoint operator in S}^ . Through a unitary transformation it 
may be represented as a self-adjoint operator in the Hilbert space ^ := L^(]R^, C'^) 
of two-spinors. More precisely, using the notation p := |p|, u)p := p/p we set 



(3) Eip) :^V^TT, Up) ^| ^^^\e{p)^^'^ ' ^ ^ ^' 
and introduce the following bounded operators on ij, 

(4) $o(p) 0o(p), $i(p) := Mp) <^ ' '^p- 

The operator $c : -^3 ^ V' ^ (*i'o(p/c) "01 'I'iIp/c) embeds ^ unitarily into 
[5]. Therefore, the form £^ defines the (two-spinor) Brown- Ravenhall operator 
in i3, 

(5) Bc[^/|x|] {ca-p + (3 - - Z/|x|) = c^E^pjc) - -Z^,(Z/|x|), 

where Uc{A) := $o(p/c) A $o(p/c) + <f>i(p/c) A $i(p/c). In the case c = 1 we 
denote this operator by Bz- Further properties properties of Bz and its relation 
to the corresponding Chandrasekhar operator and Schrodinger operator 

(6) Cz:=(p' + l)^/'-l-Z/|x|, 5z:-ip2-Z/|x| 

all realized in io, can be found in Sections [2] and [3] below and in Appendix [Cl 

1.2. Main result. We are interested in the ground state energy 
E^iZ) :=inf{£|(7/^)|VeO|, ||0|| - 1} 

of the energy form ([T]) for large atomic number Z and large velocity of light c 
satisfying Note that we picked N = Z. It was shown in [S], that similarly 
to the Chandrasekhar case [JS], the leading behavior of Ef{Z) is not affected 
by relativistic effects and, as in the Schrodinger case [39], given by the minimal 
Thomas-Fermi energy 

(7) Stf(^) :=inf{5TF(p)|pei5/3(R3), p > 0, I?(p,p)<oo}. 
The latter is defined in terms of the Thomas-Fermi energy functional 



£tf(/o) := 

where, in our units, 7tf = (37r^)^/'^/2 and 



^7TF p(x)"/3 - ^p(x) 



2 7r3 7k3 |x-y| 



is the Coulomb scalar product. By scaling, one finds Ety{Z) = Ety{1) Z'^/^. 

This paper concerns the correction to the leading behavior. For the formula- 
tion of the main result, we abbreviate the negative part of an operator by A_ := 
~^X(-oo,o)(^) and introduce for < k < the spectral shift 

(8) s{n):=n-^iY^[{B,)_-{S^)_]. 
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(We use the term "spectral shift" for s for convenience although it is used in shghtly 
different meaning otherwise.) It describes the shift of the Brown- Ravenhall bound 
state energies compared to those of the Schrodinger operator. In Section [3] we show 
that s is well-defined and discuss some of its properties. In particular, we prove that 
the function s is continuous and non-negative on the interval (0, k^] and satisfies 

(9) s(k) = 0{k^) as k ^ 0. 

We are now ready to state our main result. 

Theorem 1.1 (Scott correction). There exists a constant C > such that for 
all c > Z/k^ and for all Z > 1 one has 

(10) \e!^{Z) - Etf{Z) - (i - siZ/c)) Z^\ < CZ^^^^\ 

Put differently, Theorem 1 1 . 1 1 asserts that in the limit Z ^ oo we have uniformly 
in the quotient k = Z/c G (0, k^] 

(11) E^{Z) = EMZ) + (i - s{n)) Z^ + o{Z^). 

(We do not claim that the error Z'^'^/'^^ in (fTO|l is sharp, so we only write o{Z^) 
here.) The second term {\ — s{k)) Z^ in (fTTj) is the so-called Scott correction in 
the Brown- Ravenhall model. It does not exceed the Scott correction Z'^/2 in the 
non-relativistic model [50]. Indeed, if k = Z/c stays away from zero then there is a 
relativistic lowering of the ground state energy at order Z^ . On the other hand, in 
the non-relativistic limit c — > oo with k = Z/c — *■ 0, one recovers - non-surprisingly 
- the value of the Schrodinger case. In this case © implies 

(12) E^{Z) = Etf{Z) + ^Z^ + 0{c-^Z^ + Z^^/^^). 

The Scott correction in the Brown-Ravenhall model, however, exceeds the Scott 
correction predicted by the naive Chandrasekhar model treated in |27| and an- 
nounced in t59j . This follows from the fact that sums of bound state energies of the 
atomic Chandrasekhar operator are dominated by those of the Brown-Ravenhall 
operator, cf. the proof of Theorem 13.11 below. 

1.3. Outline of the paper. The central strategy of our paper is to compare the 
ground state energy of the Brown-Ravenhall operator with that of the Schrodinger 
operator. The latter is known up to the required accuracy o{Z'^) and the leading 
contribution agrees with the Brown-Ravenhall energy. The subtraction of the cor- 
responding ground state energies results in a renormalized effective model which 
accurately describes the energy differences and is amenable to analysis. The germ of 
this idea has been presented in the simpler context of the Chandrasekhar model [57] . 
The full blown renormalization required is developed in this paper. A virtue of our 
approach is that it leads to an explicit formula for the spectral shift which can be 
evaluated numerically. We believe it would be interesting to compare this formula 
with experimental data. 

We show that the difference between the Brown-Ravenhall and Schrodinger 
ground state energies on the multi-particle level coincides, up to the required accu- 
racy, with a spectral shift on the one-particle level. A crucial step in our analysis is 
therefore a bound on the corresponding spectral shift for rather general spherically 
symmetric potentials. This is presented in Section [31 where we show that sums of 
differences of Brown-Ravenhall and Schrodinger eigenvalues decay rather rapidly 
as the angular momentum increases. 
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In Section [2] we address various aspects of hydrogenic Brown- Ravenhall opera- 
tors. An essential feature and source of difficulties, which does not occur in the 
naive Chandrasekhar model, is the non-locality of the potential energy. In partic- 
ular, instead of the usual Coulomb potential |x|^-'^ we face the 'twisted' non-local 
operator Wc(|x|^^). Estimating the difference between the corresponding potential 
energies is the topic in Subsection 12.31 Since, in contrast to the Schrodinger case, 
the eigenvalues of the hydrogenic Brown-Ravenhall operator are not known explic- 
itly, we prove upper and lower bounds in Subsection 12.11 Our bounds are sharp 
with respect to their dependence on the quantum numbers n and I. An upper 
bound is given by the Dirac eigenvalues, a consequence of the mini- max principle 
for eigenvalues in the gap. For the lower bound we overcome the non-locality of the 
potential by a non-trivial comparison argument with a super-critical Chandrasekhar 
operator. In Subsection 12.21 we prove a new Sobolev-type inequality, from which 
we derive estimates on the eigenfunctions of the hydrogenic Brown-Ravenhall op- 
erator. The technical challenge here is to prove such a result up to and including 
the critical coupling constant. 

Finally, we present the proof of our main result, Theorem ll.il in Section HI 
For the readers' convenience we collect various facts in the appendices. Appen- 
dix [^recalls the partial wave decomposition of the Hilbert space of two-spinors. 
Appendix [Bl establishes some useful properties of the twisting operators, and Ap- 
pendix [C] collects basic facts on hydrogenic Brown-Ravenhall and Chandrasekhar 
operators. Appendix |D] fills in some details in the proof of Theorem 12.21 and, 
eventually. Appendix [E] defines the one-particle density matrix giving the main 
contribution of the energy. 



2. The hydrogenic Brown-Ravenhall operator 

In this section we set c = 1 and investigate the Brown-Ravenhall operator with 
Coulomb potential 



(13) = ^/p^ - 1 - kU{\^\-') 

in the Hilbert space = i^(K'^,C^) of two-spinors, where we recall that 

(14) z^dxp^) = $o(p)ixri<i>o(p) + $i(p)ixri$i(p) 

with <I>^ defined in (g]). In Subsection l2.1l we prove sharp upper and lower bounds on 
the eigenvalues of B,^. In Subsection l2 . 21 we prove LP estimates on the eigenfunctions 
of this operator. Technically, this is expressed as a Sobolev-type inequality for the 
massless version of which is a non- negative operator. Finally, in Subsection 12.31 
we compare the potential energy of the operator B^, namely (^ip,U{\x\^^)ip'j, with 
the corresponding local potential energy {^tp, |x|^"'^V)- Foi' comparison purpose also 
the corresponding Chandrasekhar and Schrodinger operator and occur (see 

(ED). 

According to 14 and 35j the operators B,^ and are well-defined for all k < 
with =ff = B,C and 

see also Appendix [Cl Of course, for the Schrodinger operator no upper bound on 
K is needed. 
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2.1. Estimates on eigenvalues of the hydrogen atom. In contrast to the 
Schrodinger or Dirac models, the eigenvalues of and Ck are not known explicitly. 
In order to obtain upper and lower bounds on these eigenvalues, we use that the 
spectra of -B^, and may be classified in terms of angular momenta. 

As usual write L := x x p for the operators of orbital angular momentum and 
J :~ L + icr for the operators of total angular momentum. The four operators i?^, 
J'^, J3, commute pairwise, and this also holds, if or Sk, replace -B„. This 
allows a decomposition of the Hilbert space ^ into orthogonal subspaces which 
reduce such a quadruple of operators, i.e., 

3 

(16) -^=0 %^ %i:=0%/,™. 

Here Sjj,i,m is the maximal joint eigenspace of with eigenvalues j{j + 1), of 
with eigenvalue l{l + 1), and J3 with eigenvalue m. More details concerning the 
partial wave decomposition (I16p can be found in Appendix \X[ 

We denote by bjj{K), ci{k), and s/(k) the reduced operators corresponding to 
fixed angular momenta j and Z, where, strictly speaking, we consider bj^i{K) and 
C;(k) in momentum space whereas s/(k) in position space. We refer to Appendix [Cl 
for precise definitions and further discussion. 

The main result of this subsection is that for large quantum numbers n, j, and 
I, the eigenvalues of bj^ii^K) and behave similarly to the explicitly known ones 
of the Schrodinger operator s;(k). 

Theorem 2.1 (Energies of Brown- Ravenhall hydrogen). There is a constant 
C < 00 such that for all j E Nq + \. and I — j zt ^, n G N and k G (0, k^] one has 

2 2 

(17) - C . < A„ (6,- < - 



{n + iy - " - 2{n + iy 

Here and below, we denote by Xi{A) < A2(A) < . . . the eigenvalues, repeated 
according to multiplicities, below the bottom of the essential spectrum of the self- 
adjoint, lower semi-bounded operator A. Note that —n'^{2{n + l)"^)^^ = Xn{si{K)) 
on the left hand side of (|17p is the n-th eigenvalue of the Schrodinger operator 
corresponding to angular momentum /. In particular, we conclude from (jl7[) that 
for all /i > 

(18) < tr,- ; {[B, + - [S, + < C50. 

In the proof of Theorem 12.11 we use heavily the corresponding result for the 
Chandrasekhar case, which we state next. 

Theorem 2.2 (Energies of Chandrasekhar hydrogen). There is constant C < 
00 such that for all I E Nq, n E N and n E (0, k^'] one has 

2 2 

(19) - C-^^ < A„(q(k)) < 



We break the proofs of Theorems 12.11 and 12.21 into three parts, corresponding 
to the upper bound and the lower bound for subcritical and, respectively, critical 
values of the coupling constant. 
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2.1.1. Upper hound on hydrogen eigenvalues. We begin with the Chandrasekhar 
case. 



Proof of Theorem \2.SX Upper hound. The second inequahty in (jl9[) is an immedi- 
ate consequence of the inequahty yp^ + l — 1 < P^/S and the known form of the 
Schrodinger eigenvalues in the subspace corresponding to fixed angular momen- 
tum /. □ 

Next, we turn to the Brown-Ravenhall case. 

Proof of Theorem \2.1\ Upper hound. We first recall some facts about the eigenval- 
ues of the hydrogenic Dirac operator := a ■ ■p -\- (3 — k|x|~^; see Darwin [6], 
Gordon [28] and also Bethe and Salpeter [3] for a textbook presentation. The 
following subspaces of ^^(R'^, C'*), 



1 V"'' 9{r)^j,2j-l,m{'^^)J 



reduce the Dirac operator with k e (0, 1). Under the natural identification of 
^j,i,m with L2(M+,C2) the part of D{k) in Sjj^i_m is unitarily equivalent to 



dj.i (k) 



2 _ _d_ 

r dr 
(j-0(23 + l) 




The non-decreasing sequence A„(c?j_/(k)) of eigenvalues of dj^i{K) in the gap (—1, 1) 
is independent of I and given explicitly by 



1/2 



(20) A„(d,-i(K)) 



^ (n - 1 + v/(j + 1/2)2 _ ^2^ 4.^2 



n £ 



The Dirac eigenvalues reduced by the rest energy are bounded from above by the 
Schrodinger eigenvalues: for all n, I, j, and k g (0, 1) 

(21) 1 - KidjA'^)) > - -Kisii^)). 

To show (j2ip . we use A/(j + 1/2)2 „ ^2 < + 1)2 _ ^2 < ^(n + /)2 -K2 + l_n 
and expand the outer square root in ()20p up to first order which gives an upper 
bound. 

Hence the assertion will follow, if we can show that 

(22) A„(6j,/(k)) < -1 + A„(d,, 

To prove this, we fix {j, I) and abbreviate :— X[i.oo){dj^i{0)) and A_ := 1 — A^. It 
follows from the definition of the Brown-Ravenhall operator that bj^i(K) is unitarily 
equivalent to the operator A^{djj — 1)A+ in the Hilbert space A_|_(L2(K_|_^ c2))^ 
The variational principle for eigenvalues in gaps by Griesemer et al. |30[ 131] under 
the weakened hypotheses of Dolbeault et al. [H] states that 

X„{d,AK))^ inf sup/^^^^^^^^: 0^/eyeA_(L2(M+,C2)) 

VCA+(i2(K_,.,C2)), ||/||- 
dim V = n 

Since the supremum decreases when restricted to 7^ / G V", one obtains ([22)1 . □ 
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2.1.2. Lower bounds on hydrogen eigenvalues. Subcritical case. 



Proof of Theorem \2.SX Subcritical case. Since we will reduce the Brown- Ravenhall 
case in Theorem l2.1l to the Chandrasekhar case, we actually prove a slightly stronger 
statement. As explained in ((85)) . the operators q(k) are lower bounded for alH > 1 
up to Kp > K^. 

We assume that either I > 1 and < k < or else that I = and < k < 
/ Hi . For any < (5 < 1 there exist Ms > and eg > such that 



V?TT- 1 > 



(l-(5)p i{p>Ms 
C8V^I2 \i-p<Ms. 



Denoting by Xi the characteristic function of the centered ball in with radius M^, 
and putting Xo '■= 1 ^ Xi, the Schwarz inequality implies the operator inequality 

< (1 + 6-')x^{p)\^\-'x^{p) + (1 + S)Xoip)\^\-'Xoip), 

and hence 

(23) Vp^ + 1 - 1 - A^lxj-i >x.(p) (cipV2 - (1 + S-')K\^r') xAp) 

+ Xo{p){{1~S)\p\-{1 + S)k\^\-')Xo{p). 

Now choose S as the the unique solution of the equation {1 + 5)/{l — S) — k'^ /k^ 
in the interval (0,1). Then the restrictions on k imply that (1 + 6)k < (1 — 
S)ki < (1 - S)Kf for Z > 1 and (1 + S)k < (1 - S)k'^ for 1 = 0. In any case, the 
second operator in the above sum is non-negative. The variational principle hence 
implies that the n-th eigenvalue of ci{k) is greater or equal to the n-th eigenvalue 
of Xi{p) (c5P^/2 — (1 -I- 5~^)k|x|~^) Xi(p)- Again by the variational principle, the 
latter is greater or equal to the n-th eigenvalue of csp^/2 — (1 + (5^^)k|x|^^, which 
is — const K^(n -I- ^)~^. □ 

Proof of Theorem \2.1\ Subcritical case. We assume that either j > 3/2 and < 
K < oi else that j = 1/2 and < k < k^k^ /k^. We claim that 

(24) A„(q(k)) = A2„-i(q(k) lc2) < A2„-i(6jj(k)) . 

Once we have proved this, the assertion follows easily from what we have shown in 
the proof of Theorem 12.21 above. 

To establish ([24|) we use the same notation as in the proof of the upper bound 
in Theorem 12. II By the variational principle, 

A„(6j,;(«)) 

sup inf{(/,(d,,(K)-l)/)| 11/11 = l,/eA+(i2(R+,C2)),/±/4 

A+(L=(R+,C^)) 

= sup inf{(.F,/,Q(«).Fz/) : 11/11 = l,/eA+(L2(K+,C2)),/±/4 
/i,...,/„-ie 



with J-i the Fourier-Bessel transform, see (|72p . The infimum does not increase if 
the condition / S A+(L2(M+, C^)) is relaxed to / e L'^{R+,C'^). This gives the 
eigenvalues of the operator c/(k) (8) 1^2, proving ([24]). □ 
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2.1.3. Lower hounds on hydrogen eigenvalues. Critical case. 

Proof of Theorem \2.2[ Critical case. It remains to prove that 

An(co(K)) > —const K^n^'^ 

for k^k'-^ < k < k'^. We may assume that k = and will prove that for all 
T > 

(25) 7V(-t,co(k^)) := trx(-oo,-r)(co(K^)) < const t~^^\ 

Let + Xo — 1 be a smooth radial quadratic partition of unity with Xi supported 
in the unit ball and Xo supported outside the ball of radius 1/2 about the origin. 
It was shown in [27l Eq. (19)] that the localization error can be estimated by a 
bounded exponentially decaying potential v{r) < const e^^, i.e., 

+ Xo {V^Tl - 1 - «^|x|-i - i;(|x|)) Xo- 

By the variational principle it suffices to consider the eigenvalue counting function 
corresponding to the interior and exterior term separately. The interior term is 
further estimated according to 

X^ {Vp^ + 1 - 1 - '^^Ixp' - «(|x|)) X^ > X: (|p| " ^^Ixp' " COUSt ) X^ ■ 

As shown by Lieb and Yau [41^ and explained in Corollary ID.H the number of 
negative eigenvalues of the latter operator acting in the subspace corresponding to 
Z = is finite, i.e., for all t > 

(26) Ni=o (-T, x^ (IpI - K^l^r^ - const ) x^) < const . 
For the exterior problem, we note that by the variational principle 

(27) Ni=o (-r, xo {V^Tl - 1 - K^lxj-i - z;(|x|)) xo) 

< A^/=o ("T, y/^TT - 1 - x(x)(^^|xri - v{M))) 

where x denotes the charateristic function of the support of Xo- With the singularity 
gone, the result follows as in the subcritical case. Namely, similarly as in we 
cut in momentum space according to small and large momenta. Again, by the 
variational principle, the right-hand side of ([27]) is then bounded from above by 

A^'/^ol-const T, IpI - u;(|x|)) + iV;=o(-const T,p^ - w(|x|)), 

where w(r) — const x{'>'){k^ "i"^^ + ^{r)). The first term is estimated with the help 
of Daubechies' inequality [7] 

/•OO 

Ni=o{^T, IpI - w{\k\)) < trz=o(|p| ~ < const t''/^ / w(rf/^ dr 

Jo 

with the latter integral being finite. For the second term we estimate w{r) < 
const r~^ and use that 

iV;=o(-T,p^ - const |x|"^) < const r^^^^. 

This concludes the proof of Theorem 12.21 □ 
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Our proof of Theorem 12.11 in the critical Brown- Ravenhall case is based on a re- 
duction to the Chandrasekhar case. The next lemma compares the number of eigen- 
values of the critical operators 61/2, with those of the two operators ci'{k^) 
with /' = 0, 1 and critical coupling constants = 2/n and = tt/2, cf. 



Lemma 2.3. There exists a constant such that for 1 — 0,1 and all t > one has 
N {~T, 6i/2,i(K^)) < const [N (-T, co(k^)) + N {-r, ci(«;f ))] . 

Proof. We start with the observation that + (nf)^^ — 2{k^)^^. Using the 

explicit form of the reduced operators (cf. Appendix [Cjl . this implies the identities 



&l/2,o(«:-^) = [{1^0) Vofoo,0'/>0 + (k? )~Vl&l,l</'l 

(28) / ~ ~ X 

^i/2,i(«:'^) = «^ (^(k? )"Vi6o,i'/'i + )"Vo&i,o0oj , 

where the operators bi^i, are defined in L^(M+) through quadratic forms 



{f,bu.f)--^ ^^^|/(p)|^dp-«;^y^ f{p)k'r{p,q)f{q)dpAq. 

In case v = 1 it hence follows from 2(l)i{pf < 1 that (/,6;,i/) > {f,ci{Kf)f). In 
case = we use the inequality 



(29) {E{p) ~ l)^o(p)-' > Vp'+4 - 2 - 2(i?(p/2) - 1) 

which is most easily seen by writing both sides in terms of E{p). It implies 

(/,&i,o/) >2(w/,Q(«;p)«/) 

where the unitary scaling transformation u is defined through (uf){p) := %/2/(2p). 
The proof is completed by the variational principle. □ 

We are now ready to give a 



Proof of Theorem \2.1\ Critical case. The previous lemma implies that it suffices to 
show that for I = 0, 1 

N (-T, q(kP)) < const T"^/^. 

In case / — this was established in (I25p . and the case I — 1 follows similarly with 
the analogue of (pS)) given in Corollarv lD.il □ 

2.2. Sobolev inequality for the critical Brown-Ravenhall operator. Having 
studied the eigenvalues of -B^ in the previous subsection, wc now turn to integra- 
bility properties of its eigenfunctions. The L'-norm of two-spinors is given by 

1/9 

1^11,:= ( / |^(x)rdx; 



where the modulus, | • |, refers to the Euclidean norm in C^. For q = 2 wc drop the 
subscript. We aim at proving the following 

Theorem 2.4 (L'^-properties of eigenfunctions). Let 2 < q < 3. There exists a 
constant Cq < 00 such that for any k E (0, k^] and all ip G 13(5^) with (^/;, B^V') ^ 
one has ijj G L"^ with 

(30) \m,<Cq\m 
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Note that (|30|) applies, in particular, to eigcnfunctions of B^, corresponding to 
negative eigenvalues. The proof of Theorem 12.41 which is spelled out below, relies 
on a Sobolev inequality for the massless atomic Brown-Ravenhall operator in ^) 
given by 

Bl") := \p\-^{\^\-'+u;^-a\^r ■ ^) ■ 
This operator is bounded below (in fact, non-negative) if and only if k < . 

Theorem 2.5 (Sobolev inequality). For any 2 < q < 3 there exists a constant 
Cq>0 such that for all G Q(B^b), 

(31) Ml<C, \m'^'-'\ ^ = 6(i-|). 

It is illustrative to compare pTjl with the 'standard' Sobolev-Gagliardo-Nirenberg 
inequalities, 

(32) ll^ll^ < C; i^P, M'^'~'\ ^ - 6(i - i), 2 < g < 3, 

see, e.g., [38l Thm. 8.4]. Hence Theorem 12.51 savs that, if the endpoint g = 3 is 
avoided, an inequality of the same form remains true after subtracting the maximal 
possible multiple of Z/^(|x|^^) from |p|. Moreover, one can show that ([3T|l does not 
hold with (7 = 3, not even if the L'^-norm is replaced by the weak L'^-norm. 

Note that if k < then ^ with 5^°^ instead of Bf^ follows from §^ - 
but with a constant that deteriorates . The main point is to derive an 

inequality which holds uniformly in k up to and including the critical constant. Our 
proof is based on the somewhat surprising fact that the Brown-Ravenhall opera- 
tor with coupling constant can be bounded from below by the Chandrasckhar 
operator with smaller coupling constant k*^. 

Before we start the proof of ([3T|) . we provide the 

Proof of Theorem \2.4\ The Sobolev inequality (|3ip implies 

<C,\\B^^'>-B4'Uf. 

Tlx showed [Ml Thm. 1] (see also Balinsky and Evans \T) that the difference 
_b1°^ — _Bk extends to a bounded operator with norm uniformly bounded for any 

Ke(0,K^]. □ 

2.2.1. Comparison of critical operators. The first step in the proof of the Sobolev 
inequality (PT|) is a comparison of B^^ with the massless atomic Chandrasckhar 
operator in i^, which is given by 

:= \p\-n\^\-\ 

It is bounded below if and only if k < k*^. As discussed in Appendix [Cl those parts 
of B^i'^ and ci^^ in the subspace S^j,i,m, are unitarily equivalent to operators fe^""* (k) 

and cf^\K) in L^(Ili+), which depend only on j in the Brown-Ravenhall case and 
only on I in the Chandrasckhar case. For the comparison argument it is important 
to note that the reduced operators b^^\K) and c[°^(k) are lower bounded for k up 
to and including the critical coupling constants and Kp respectively. They are 
defined in ([M]) and, as is explained there, exceed and k*^, if j > 3/2 or Z > 1. 
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We begin by observing that all the critical operators b^^\K^) and c|"''(Kp) have 
the same 'generalized ground state', namely p. The corresponding ground state 
representation formula (in momentum space) is given in 

Lemma 2.6 (Ground state representation). // / G £l{b''p [k^)) and g{p) = 
Pf{p), then 

(33) if,bf\.f)f)^^ r r\gip)-giq)\^kfim + ^))^^. 



Similarly, if f E 0(c|'''' (Kp)) and g(jp) —pf(jp), then 



^ 1-^ 1-^ dpdq 



(34) (f,cf\^r)f) = ^ / l5(p)-5(g)pA:F(i(2 + l)) 

Jo Jo PI 



where and ki are defined in ((81 



Proof. We write k for one of the functions kf or fcp and n for the corresponding 
constant or Kp. Expanding the square, we find 

1 f f ./„M2,.l/p , a^^dpd(7 



, . |5(p)-5(g)|^fc(i(H + i)) 



oo 

l5b)P 



'p|/b)P (/ + f )) dp - f(p)K\{l + §))/('?) dpdg. 



By definitions ([84)1 and ([85]) of k we have 

which implies the assertion. □ 

Now we bound B^^^ from below by C^^c ■ 

Lemma 2.7 (Comparison of critical operators). There is a positive constant 
such that for any ip G 0(5^^^) (H ^^^2 i 

(35) ( V, Bf, ^) > const ( V, C^c' ^) 

An inequality of the form ([SS]) cannot hold in the subspace ^1/2,1, since the right 
hand side is bounded from below by a constant times (V', |p|?A) while the left hand 
side is not. 

Proof. By orthogonality it suffices to prove the inequality on each subspace 
First let {j,l) = (1/2,0). We may also fix m = ±1/2 and choose tp G ■fti/2,o,m- Its 
Fourier transform is of the form ■0(p) — P^^/(p)^i q ('^p); see Appendix El By 
the massless analog of ([55| one has 

(V,Bi°j^) = (/,fe[%(«^)/). 
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Setting f{p) —: pg{p) we obtain in view of Lemma 12.61 

^ Jo Jo PI 

> (1 + rr \9{P) ~ 3(9)Pfco^(i(? + f )) 

^ Jo Jo PI 

> (l + (2/7r)2)-i(/,4°)(^^)/) (l + (2/7r)2)"i(V>,cl°jV) 

Here we used that > Qi and the massless analog of ([5^ . This proves the assertion 
on the subspace i3i/2.o- Now assume that i/j e {Sji/2,0 ® -^1/2.1)^ ^nd note that 

|p|>^(|xri+c.p.^|xria;p.^) 

on that space. Here we used that is monotone increasing in j, see Appendix ICl 
We conclude that 

^3/2 '^3/2 

proving the assertion. □ 
2.2.2. Proof of the Sobolev inequality. We are now ready to give a 
Proof of Theorem \2.5\ By scaling, pip is equivalent to the inequality 

\mi<c', ((v,si°i^) + iiv^ip). 

This, together with the triangle inequality, shows that it is enough to prove the 
inequality separately on the subspaces ^1/2,1 and ^^^2 v latter subspace, 

the claim follows immediately from Lemma 12.71 above and the Sobolev inequality 
for the critical Chandrasekhar operator 26, Corollary 2.5]. We now reduce the 
claim for the subspace ^1/2,1 to that for f)i/2,o- For this purpose, we note that the 
helicity operator H = • cr, cf. (|73p. commutes with b]^^ and, by ([7^ . maps 
into S)j,2j-i. Hence if G -^1/2,1 then by the Sobolev inequality on fii/2,0 



\il^,B'^>i:J + llVi' = (^M^,BI'>M^'J + IIHV^II^ > const HHV^H^. 
By Lemma |BT] the hehcity is bounded on (R^ , ) . □ 



2.3. Estimates on the electric potential. The goal of this subsection is to 
compare twisted and untwisted electric potentials. We begin with an estimates for 
point charges and then turn to smeared out charges. 



Lemma 2.8. Let I > 1 and ijj e Then 



const 



(36) KV^dxI-i-Z^dxI-i))^)! <-^(V,p^^). 

Proof of Lemma \2M By orthogonality it suffices to prove the assertion for -0 S 
^j,i,m- Its Fourier transform is of the form ^/'(p) = /(p)p~^r2jy_m(wp), cf. Appendix 
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lAl and we compute similarly as in 

(^^,(1x1-1-^^(1x1-1))^) 

= i ^pW)j^ M /(g) {[1 - Mp)Mq)] Qi + i 

-0i(p)0i(g)g2,-z(i(§ + f))} 



with 



Jo Jo „_n ^ ^ 



dpfip) / dg/(g)0i(p)^i(g) 
Jo 

We estimate these terms separately. For the first term we use (fTf]) and ((78)) together 
with Abel's argument to turn Hermitian integral operators into multiplication op- 
erators by means of the Schwarz inequality (see also [HI Ineq. (6.9)]). Since the 
Qi are positive, we obtain 

- r w r ill ' ^„ ^^^^^^ ~ ^ ') 



<- 



dp 



\fip)\' 



2 /-oo 



E{p)^ Jo 



dp 



E{p)^^ 



dq 



dq 
9^ 



ip~<ifQi{\{i+i)) 



{l-qfQi{\{q + q-')) 



We now use the bounds p^ /E{pY < P^ and, for g > 1, (1 — q)^ < 9^ — 1 which yield 



dq 
q" 



{^-qfQiiW^q-^)) 



< 4 



dq 

1 q 

dx Qi{x) 



.{l^q^QliHl + l^')) 
4 



where the last step involved p^, 324(18)]. Thus, 



2l{l + l) 



dppV{p)\' = 



21(1 + 1 



■(^,p'V')- 



We estimate the term A2 similarly by the Schwarz inequality, 



A,<2 



<4 



dp\f{p)\'\Mp)? 



dq 



P 



dq 



dp\fip)\V I j\Qi{Hi + i-'))-Q^^-dHi + i-'))\ 



Due to the pointwise monotonicity ([57)1 the difference inside the modulus is of 
definite sign. Without loss of generality, we may therefore assume 2j = 2/ + 1. 
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Using the integral representation (jC.ip we can bound 



da; 1 



y/x'^ - 1 Vl - 2xz + z2 

Adding the estimates for Ai and A2 we arrive at ((36|l . □ 

Note that our proof shows that one can choose different powers of |p| on the 
right hand side of ((36| . 



Lemma 2.9. There exists a constant such that for any electric potential v of a 
spherically symmetric non-negative charge density 

\{ip,{v ~U{v))tJj) \ < const w(0) (V',P^'/')- 

Proof. We denote by r : K."^ — s- [0, 00) the sphericaUy symmetric, non-negative 
charge density corresponding to v, i.e., f(x) ~ /'''(x — y) |y|^^dy. The Fourier 
transform of r obeys the estimates 



l^(P)l = \/— 2/ ''|sin(|p|r)r(r)|dr < 



(2^)3/2 

By Fourier transform the scalar product on the left side of the assertion becomes 
(^, {v - U{v)) II V'(p)* jp^__~jJ (1 - <fo(p)<i>o(q) - <i>i(p)$i(q)) ^(q)dpdq. 



Using Lemma IB .21 below we estimate the absolute value of the preceding expression 
from above by two terms, Bi and B2. The first term can be further bounded as 
follows. 



Bi= const J J |f(p-q)||^(p)||V'(q)|dpdq 

< const t;(0) I I (^M^'^ \^\'^'^ 

< const i;(0) y"|-0(p)|^p^dp, 

where we use the Schwarz inequality in the second step. The second term is esti- 
mated similarly 



52 = const // |f(p~q)|^^^|^(p)||V;(q)|dpdq 



< const .(0) /dp|^(p)l7(^) P^dq 

< const i;(0) yiV'(p)P|ppdp. 



□ 
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3. Spectral shift from Schrodinger to Brown-Ravenhall operators 

The main theme of this section is the (integrated) spectral shift, i.e., the differ- 
ence of sums of eigenvalues of the Brown-Ravenhall or Chandrasekhar operator 

B[v] -.^ y/p^ + l-l -U{v), C[v]:^ E{p)-l-v, 

(cf. (fT4|) ) and the Schrodinger operator S[v\ := ip^ — v, all acting in the Hilbcrt 
space Sj of two-spinors. We have set c = 1. 

Concerning the potential w : M we will always assume that the above oper- 

ators can be defined through the Friedrichs extension starting from ©(M'^, C^). For 
example, the condition < w(x) < k'^ |x|~^ with ^ ~ B, C (cf. (HH)) ensures that 
the Brown-Ravenhall, respectively the Chandrasekhar operator are well-defined and 
bounded from below (see [14] and [35]). 

We assume throughout that the potential v is radially symmetric which allows 
us to investigate the spectral shift on each subspace S)j,i in the decomposition ([16]) 
separately. We write Ajj for the orthogonal projection onto Sjjj. For the reduced 
traces we use the notations 

tTj.i{A) := tr(Aj- zA), trj{A) tr^j+i/2(A) + trjj„i/2(^). 



3.1. Estimate on the spectral shift. One of the key observations in our proof 
of the Scott correction is that the spectral shift between the one-particle Brown- 
Ravenhall and the Schrodinger operator decreases sufficiently fast for high angular 
momenta. 

Theorem 3.1 (Spectral shift: Brown-Ravenhall case). There exists a con- 
stant C < oo such that for any k < , any v : [0, cxj) — > [0, cxd) satisfying 

(37) v{r)<Kr-\ 
any ^ > and any j G Nq + 1/2 one has 

(38) tr, {[B[v]+^i]_ - {S[v]+^i]_) < C j-\ 

We derive this result from a corresponding theorem for the Chandrasekhar oper- 
ator. For a proof of the latter we need to strengthen f27!, Thm. 2.1]. In particular, 
we need to consider C\v\ for potentials v satisfying ([37]) also in case < k < . 
Those operators are not densely defined in the Hilbert space i3. However, according 
to ([85]) below, they are densely defined in the subspaces Sjj^i with j > 3/2. Another 
new aspect is that we trace the dependence on the coupling constant. 

Theorem 3.2 (Spectral shift: Chandrasekhar case). There exists a constant 
C < CO such that for all I E Nq, j ~ I zt ^, for all k satisfying 

r zfl = 0, 

-{ ^fl>l, 

for all fi > and for all v : [0,oo) — > [0,cx)) satisfying ()37|) . one has 



(39) < tY,,i {[C[v] + - [S[v] + /i]_) < C 



A 



(1 + ^) 
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One of the key points to be appreciated in the above theorems is an effective can- 
cellation in the differences in (|39|) and ([38]) . This can already be seen for Coulomb 
potentials v(r) — nr"^, where 

2 ^ -1 



2 ^ (n + O^' 

n— 1 ^ ^ 

which does not decay at all as j ^ oo. Moreover, for fixed j and I the above trace 
vanishes only like 0. It is rather remarkable that such cancellations occur 

uniformly for all attractive potential v satisfying (j37p . 

The following proof of Theorem 13.21 follows the ideas of [27j Thm. 2] . It is not 
only included to render the paper self-contained, but also to establish the above 
mentioned improvement, which arc important for the present paper. 

Proof of Theorem \3.Sl We note that both traces tr^ ^ [C[v] + and tr^ ^ [S[v] + 
are finite. This follows by the variational principle from the case v{r) = Kr~^, cf. 
Theorem 12.21 in the Chandrasekhar case. Thus, for / < 3 say, it is enough to show 
the claim for k in a neighborhood of 0. More precisely, we can assume k < -^(1 + ^) 

which covers all K< for l> 3. 

Moreover, by an approximation argument it is sufficient to consider /i > and 
bounded potentials v, cf. p7] . 

We denote by djj the orthogonal projection onto the eigenspace of C[v] cor- 
responding to angular momenta j, I and eigenvalues less or equal than — /i. The 
identity 

(40) ^ ^ 1^2 

and the variational principle (cf. 38, Thm. 12.1]) imply 

(41) < 2tTj,i {[C[v] + /i]_ - [S[v] + fi]_) < tr [C^dj,i] . 



Using the eigenvalue equation and the bound (j37p on the potential we estimate this 
term further as follows. 

tr [C^d.^i] < trj-z [C[v]f_ + tr [v^d.^i] < tr.j [C^t + kHt [\x\-^d,^i] . 

Using Hardy's inequality and (I40p 

tr [\^r^dj,i] <{l + i)"'tr [p^d,,i] = {l + i)"' (tr [C^dj,i] + 2tr[Cod,.i]) . 

Since k < I + ^, the last two estimates may be summarized as 

(42) tr [C',d,,] <(^1- (7^) (tr„; [C.]^ + -^j^ tr [Cod,, 

We shall estimate the two terms on the right hand side separately. From p71 
Lemma 3] we recall the following angular momentum barrier inequality on 

(43) Co > 2Kr-\[r<RM}, Rii^i) = j^il + i)'. 
(Here we use that k < -^{l + ^)-) This implies 

tr [Codjj] < Ktr [|x|"^djj] < ^ tr [Codj,/] + ^ tr [widj^i] 

= |tr [Cod,v]-itr [C[wi]d,,i] 
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where wi(r) :— Anr ^X{r>Ri{K.)}- Hence, using the variational principle followed by 
Daubechies' inequality [7] (cf. also [27l Prop. 1]) 

/ />00 /"OO 

tr [Codj,i] < tr j,i [C[wi]]_ < const {21 + 1) f j wi{r)^^^ ^ J "'''^''^^ 
(44) < const K^. 



In order to estimate the first term on the right hand side of (|42|) we use (|43l) to 
obtain on fyj^i 

C« > iCo - Kr"^X{r>fl,(K)} > \C[wi\. 
with wi as above. Hence again by Daubechies' inequality 

tr,, ; [C«]i < const (2/ + 1) (^j^ wi{rf''^ ^ / "''''^''^ " ^ 

Combing this with ([44]) , (|42]) , and (|41]) completes the proof. □ 

Having finished the proof of Theorem 13.21 it is easy to give the 

Proof of Theorem \3.1\ Since the trace tr,- [B[v\ + is finite according to Theo- 
rem [2lT] we may assume that either k < and j = 1/2, or else that j > 3/2. In 
this case, the claim essential boils down to Theorem 13.21 To see this, we note the 
identity 

(45) B[v] = U{C[v]) - i {U{prC[v]U{p) + U{p)C[v]U{pr) 

involving the unitary operator L/(p) :— <I>o(p)+«<I'i(p) (see also Equality (l45l) 

as well as the unitarity of U{p) are easily derived from the fact that <i>Q(p)-|-<i>^(p) = 
1. 

Even if v satisfies (|37l) only with a. kP < n < ^ identity (j45| remains valid 
on all subspaces io,- with j > 3/2. Hence by the concavity of the sum of negative 
eigenvalues [65 of B[v] + /i one has for any /-i > 

tr, [B[v] +^^L<l tr, [U*{p)C[v]U{p) + + I tr, [U{p)C[v]U* {p) + m]_ 

(46) -tr,[CH+Ai]_. 

By ([5^ the trace in ([55)1 is thus bounded from above by 

tr,- {[C[v] + fi]_ - [S[v] + fi]_) < const n^j"^, 
as claimed. □ 

3.2. Properties of the spectral shift. In this subsection we discuss some prop- 
erties of the spectral shift s{k) defined in ([5]). 

Lemma 3.3 (Properties of the spectral shift). The spectral shift s is a con- 
tinuous, non-negative function on {Q,k^] satisfying s{k) — 0{k'^) as k 10. 

Proof. According to and Theorem 13.11 one has 

< Sj(k) K^^tr, {[B^]_ - [S^]_) < const 

Therefore the sum s{k) = converges, is non-negative and satisfies the 

claimed asymptotic estimate as k | 0. By the mini-max principle each eigenvalue 
depends continuously on k. Thus the continuity of their sum follows from the 
estimates in Theorem l2.1l and the Weierstrafi criterion for uniform convergence. □ 
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4. Proof of the Scott correction 

The strategy of the proof of the main results is similar to the one used for 
the Chandrasekhar operator [27|. We employ the Schrodinger operator as a reg- 
ularization for the relativistic problem, i.e., we will use it to eliminate the main 
contribution to the energy (the Thomas-Fermi energy) and focus only on the en- 
ergy shift of the low lying states. For these the electron-electron interaction plays 
no role and the unscreened problem remains. We define 

E^{Z) :=inf{f|(V^)|V'GQ|,||V'll =1} 
to be the ground state energy in the Schrodinger case, 



£'nW ■■= ( V', 



v=l ^ ' l<fj.<v<N 

It is defined on Qff := i3f,n6(R^^, C^"), where Sjff := A^Li ^ is the Hilbert space 
of anti-symmetric two-spinors. We recall that we suppose neutrality, i.e., N = Z. 

The asymptotics of the Schrodinger ground-state energy up to Scott correction 
reads [5U] 

(47) E'^{Z)=ETFiZ) + ^Z^ + OiZ^^^^^). 

For our purpose this remainder estimate is sufficient. However, even the coefficient 
of the Z^/^-teim in the asymptotic expansion is known pUl [ril [T^ IT51 ITBl 

[mini. 

Our main result. Theorem 1 1, li will follow from (j47p if we can show that in the 
limit Z — > oo the difference of the Schrodinger and Brown- Ravenhall ground-state 
energy satisfies 

(48) E^{Z) - E^{Z) = s{Z/c) Z"^ + 0{Z^'^'^'^) 

uniformly in k = Z / c G (0, k^]. We break the proof of this assertion into an upper 
and lower bound. 

4.1. Upper bound on the energy diflference. The Thomas-Fermi functional 
([7]) has a unique minimizer qz, the Thomas- Fermi density (Lieb and Simon j39j). 
It scales as qz{^) '■= Z^gi{Z^/'^x). We set 

(49) <^tf(x) ■.= Z\^\-^- f ^^dy, 
the Thomas-Fermi potential, and 

iTF(x) := / -. rdy. 



Kx) := / 



X- 



<;-y|<-Rz(x) 

the exchange hole potential. Here i?z(x) is defined as the (unique) minimal radius 
for which i|x-y|<_R,z(x) Sz{y)dy = \. The corresponding one-particle operators - 
self-adjointly realized in - are 

S'tF — S[(j)TF + -^Tf], Stf = Bc[(i>TF + ^Tf]- 

Here we use a notation analogous to that in ([5]). 

We shall express the many-particle ground-state energies E^{Z) and E^ [Z) in 
terms of quantities involving the above one-particle operators. In the Schrodinger 
case, this was achieved in [50l[53] in terms of the Thomas- Fermi potential </>tf- Our 
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point in the proof of the following proposition is to replace (jjTF by the exchange 
hole reduced potential (jjTF + Ltf- 

Proposition 4.1. Let J := [Z^/^] + i. Then, as Z oo, 

J-l Z+l/2 

(50) £;«(Z) = - J2 tr,[S[Z\^\-']]_- tr,[STFL-D{gz,gz) + 0{Z*y^^). 

j = l/2 j=J 

Since 4>tf + Ltf has a Coulomb tail, the trace tr^ [S'tf]_ is finite for each j, 
but not summable with respect to j. It is therefore essential to restrict the second 
sum to a finite number of angular momenta. However, the value of the cut-off, 
j < Z + l/2, is not chosen optimally here, since for our argument it is largely 
arbitrary. 

Proof of Proposition \4-l\ According to the correlation inequality [42] 

z+l/2 

E''{Z)>- J2 t^j[STF]^-D{gz,gz). 
i=i/2 

Note that the Z electrons can certainly be accommodated in the first Z angular 
momentum channels (which is a very crude bound). Estimating 0tf + ^tf from 
above by the Coulomb potential for small angular momenta, we obtain 

J-l 2+1/2 

(51) E^{Z) > - ^ tr, [S[Z\k\-']]_ - ^ tr, [StfL - D{gz, gz)- 

j=l/2 j=J 

Moreover, see |50l [53] , 

,7-1 oo 

E'^iZ) <-Y. t^J- [S[Z\^r']] _ - E trj- [^[-^tf]]- - D{gz, gz) + const Z^'/^\ 

3 = 1/2 j=J 

Hence it suffices to prove that 

2+1/2 oo 

(52) - iij [Stf\_ > - E [-^[^tf]]- - const Z^l^ 

J=J 3=J 

(Note that the lower bound in [57] contains an error by estimating [571 Equation 
(43)] to generously. Really, only the first Z lowest negative eigenvalues need to 
occur on the right hand side instead of all. In particular, there will be never more 
than Z total angular momentum channels occupied. This fact is taken into account 
here yielding a suitable lower bound. The problem in [27] can be circumvented in 
exactly the same way.) We decompose Ltf = L< + L> where 

L< = X{|x|<7?}-^TF, L> = X{|x|>fi}iTF, 

with a constant R (independent of Z) to be chosen below. For e > to be specified 
later we estimate using the variational principle for sums of eigenvalues 

(53) tr, [5tf]„ 

< trj(i(l - 2e2)p2 _ </,Tp)_ + £2 trj.(ip2 _ £-^L<)_ + trj(i/ _ £-^Ly)_. 
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By the subsequent lemma the first and main term is bounded according to 

2+1/2 oo 

tr,(i(l - 2e^)p' - - J2 - -^tf)- 

< tr(i(l - 2e^)p^ - 0tf)- - tr(ip2 _ < const e^Z^^^. 

For the second term on the right side of ((53|) we use the Lieb-Thirring inequality [40] 
to obtain 

72/3 



< const e"^ J L<(x)^/2 dx < const e^^Z^ 



In the last inequality we used a bound of Siedentop and Weikard ^53^, Proof of 
Lemma 2]. It is at this point that R is chosen. The penultimate inequality in 
[53l Proof of Lemma 2] asserts after scaling that i>(x) < const |x|^^. Hence by 
comparison with the exact hydrogen solution 

Z+l/2 Z+1/2 

trj(ip2 - e"2^>)_ < ^ tr^-^p^ - e^^const |x|"i)_ 
j=J i=i/2 

2+1/2 oo 

— 2 X — ^ X — ^ J ' —9 

= const e > > -—t: < const e Z. 

, "^An + j- 1/2 2 - 

j=l/2 n=l ^ •' I ' 

Choosing e = Z^^l'^ all the error terms arc 0(Z^I'^\ proving ((52)) . □ 

In the previous proof we used 
Lemma 4.2. For all < e < 1/2, as Z — > oo, 

(54) tr(i(l - £2)p2 _ < lp2 _ ^ ^Qjjg^ £2^7/3, 

Note that there are only a finite number of eigenvalues, since <f)TF decays like 
|x|-4. 

Proof. Let c?|.p be the projection onto the negative eigenvalues of ^(1 — e^)p^ — 0tf- 
Then, by the variational principle 

(55) tr(i(l - e2)p2 _ _ ^^(1^2 _ 

< -trd^F(5(l - 0tf) +trc?^F(5P^ - 0tf) = ^tr^TEP^- 

Hence the claim will follow, if we show that tr d^pp^ < const Z^/"^ . Note that dj^p 
depends on both e and Z, and by rescaling one may get rid of the e dependence 
at the expense of changing Z. We may therefore assume that e = and write 

^TE = ^TE- 

Thus, it remains to prove 

(56) tr dxEP^ < const Z^/^ 

Note that this says that the kinetic energy is bounded by the order of the total 
energy tr dTF{^p'^ — (/)te), which is well-known to be of order Z''/'^. Using that 0te 
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is bounded by a constant times min{Z|x|^^, |x|^^} (see [39]) we get for any R> 
itrdxFP^ < trdxp'/'TF 

3/5 

< const 



U ^^^^(Z|x|-l)5/2dxj dTF(x,x)5/3dx^ 

J t^TF(x,x) dx^ 



The Cwikel-Lieb-Rozenblum inequality (for a textbook presentation, see, e.g., [55] 1 
guarantees that 

y dTF(x, x) dx < const J 0TF(x)'^^^dx = const Z. 
Moreover, by the Lieb-Thirring inequality [10] 

J (iTF(x, x)'"^/^ dx < const trdxFP^- 
We can estimate for any S > 

(\ "^/^ 3/5 
y^ ^ ^^(Z|xrl)5/2dxj (^yrfTF(x,x)5/3dx^ 

< const Zi?^/^ (trdTFP^)^^^ 

< (5trdTFP^ + const (5^3/2^5/2^1/2^ 

In summary, we have shown that 

(i - const S) trdxFP^ < const (5-3/2^5/2^1/2 _^ ^-4^^ ^ 

Choosing S small (of order one) and R ~ Z^^l^ we obtain □ 

Next, we bound the many-particle ground state energy of the Brown-Ravenhall 
operator from below by one-body quantities which match the corresponding quan- 
tities in the Schrodinger case ([St 



Lemma 4.3. For all J e No + 1/2 and Z eN 

.1-1 z+1/2 

E^{Z) >-Y. t'-J- [Bc[Z\^\-^]]_ - t'-^- [Bty]^-D{bz,Qz). 

j=i/2 i=J 

Proof. This follows by the same argument leading to (|5T]) . □ 

We are now ready to give a 

Proof of Theorem \l.l\ - first part. Choosing J — [Z^/^] -I- i and combining Propo- 
sition 14.11 and Lemma 14.31 we obtain 

(57) E^{Z)-E^{Z)< tr,([B,[Z\^\-^\]_-[S[Z\^\-^]]_) 

J = l/2 
^+1/2 

+ tr, ([i?TF]_-[^TF]_)+0(^'''/''). 
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We note that by scahng x i-^ x/c, the operators S'[Z|x|^^] and _Bc[^|x|^^] are 
unitarily equivalent to the operators Z'^kT^Sk. and Z^k^^B^ where k = Z/c. Sim- 
ilarly, 5tf and -Btf are unitarily equivalent to the operators Z^k~^S'[k|x|~^ — Xc] 
and Z'^k~'^B[k\x.\~^ — Xc] acting in 9), where 

Xc(x) c M -dy. 

J|x-y|>ci?,z(c-ix) |X-y| 

This implies that the first two terms on the right-hand side of ([57|) , which we denote 
by ^i{Z, c) and 'S2{Z, c), can be rewritten as 

Si(Z,c)=Z2k-2 ^ tr, {[B,]_ - [S,]_) , 

J = l/2 
2+1/2 

S2(Z,c)=Z2k-2 tr, ([i3[«|xri-xc]]_- [5[«|xri-xc]]_^ 

3 = J 

Inequality (jf 8p and Theorem 13 . f I guarantee that the terms in the first sum are non- 
negative and that the terms in both sums are bounded from above by a constant 
times independently of Z and c. Therefore, the first sum can be bounded 

from above by an absolutely convergent series, 

oo 

i=i/2 

By the same token 

oo 

^2{Z,c) < const Z^K^^r' = 
j=J 

uniformly in c. This concludes the proof of the upper bound on the energy differ- 
ence. □ 

4.2. Lower bound on the energy difference. Similarly to [50] we define one- 
particle density matrices and d^ on Sj as sums 

(58) d*=d*+dy, if^S,B. 

The contribution of small total angular momenta, — ^i^i^df, is defined in 
Appendix IE. II It comes from the eigenspinors of the atomic problems. The con- 
tribution of large angular momentum, (i> = X^z^l^'i defined in Appendix IE. 21 
It corresponds to the Macke orbitals of f50l and, in particular, coincides for the 
Schrodinger and Brown-Ravenhall case. The angular-momentum cut-off L will be 
chosen in a Z-dependent way, namely. 

Important properties of the density matrices, whose construction is explained in 
more detail in Appendix [El are: 

• The densities 

p#(x) := trc2 (d#(x,x)) , pf (x) := trc2 (df (x,x)) , 

KL 1>L 
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of d"^, df , and d > are all spherically symmetric. 
• The dimension of the ranges of the density matrices d^ and d^ is at most 
Z, in particular trd'^ < Z. Moreover, 

(59) tidf = J p*{x)dx^2{2l + l){K -I), 0<1<L, 

with K = [const Z^^^] and a suitable constant. 

For a lower bound on the ground state energy in the Schrodinger case, we recall 
from [5D] and [571 Proposition 4] the following 

Proposition 4.4. For large Z, 

E^{Z) = tr [S[Z\x\-^] d^] + D{p^, p^) + 

To obtain an upper bound on the ground state energy in the Brown-Ravenhall 
case, we use the reduced Hartree-Fock variational principle. It involves the density 

p^(x) :=trc2 (Z^,(d^)(x,x)) 

of the twisted density matrix Uc{d^). 
For further reference, we also set 

:= trc2 {Uc{df){x,x)) , p^,<(x) ;(x), p[/,>(x) := ^py.iW- 

KL 1>L 

Applying to ([T]) the Hartree-Fock variational principle - in the strengthened version 
of Lieb [37] (see also Bach |I]) - and omitting the manifestly negative exchange 
energy we arrive at 

Proposition 4.5. For all Z and c, 

E^{Z) < trmZ\x\-']d'']+D{p§,p§). 
Combining Propositions 14.41 and 14.51 we find 

Ef{Z)-E''{Z) 

< tr[Bc[Z\x\~^]d^] - tr[S[Z\x\-^]d^] + D{p§ ~ p^ , p§ + p^) + const Z^^/^^. 

Now we use the inequality p^ > 2c^{E(p/ c) — l) for the kinetic energy corresponding 
to dy. Morover we remark that D{pu,> — P>tPu<) — ^3 D{f,g) < 0, if 
/ < < .9. This yields 



(60) E^{Z)-E^iZ)<tr B,[^^]d^ - tr 



tr 



|x| 



D{pu,> - P>.Pu,> + P>) +2 D{p§^^,p§ + pS) +const Z''^^^ 



As we shall see, the first two terms will yield the Scott correction. In the following 
subsections we prove that TZi, TZ2, and 7?.3, are relatively small remainder terms. 
Hence, we wish to control the effects of the twisting operation Uc, which stems from 
the electronic projection, on the electrostatic Coulomb energy. 
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4.2.1. Controlling the electron projection for high angular momenta. Our task in 
this subsection is to prove that for large angular momenta, the twisted and un- 
twisted electrostatic energy are asymptotically equal. 

We start by comparing the electric potential energy with or without electron 
projection for large angular momentum. This will imply that the term TZi in (j60|) 
is relatively small. 

Lemma 4.6. In the limit Z ^ oo one has uniformly in k — Z/c G (0, k^] 
J (p>(x)-p^,>(x))^=tr[(|xri-Z^,(|xri)) d>] ^ 

Proof. Let {ipa} stand for the Macke orbitals building up rf> which we label by 
a — {j,l,m,n); see (j95p and preceding equations in Appendix IE. 21 By the scaling 
X I— > x/c one has the relation 

where il^a\'^) c^'^/^i/'q(x/c). Assuming that a corresponds to a fixed (large) 
(j, /) we may use Lemma 12.81 to estimate the right-hand side by a constant times 

Using that Z/c < we obtain the estimate 

B °° 1 

tr[(|xri-Z^,(|xri))d>] < const ^ ti;,i [p'd^] . 

l=L j=l±l/2 

The proof is completed using Lemma IE. II from Appendix IE. 31 □ 

Next, we estimate the difference of Coulomb energies corresponding to large total 
angular momenta. This shows that the term TZ2 in (|60p may be neglected. 

Lemma 4.7. In the limit Z oo, one has uniformly in k ^ Z/c G (0, k^] 

7^2 = D{pu,> - p>,p> + pu,>) = o{z''i-^). 

Proof. We define v :— (p> + pu,>) * I ' 1^^ to be the electric potential generated by 
P> + Pu,> which is obviously spherically symmetric and obeys 

z;(0) = tr [d> (1x1-1+^^,(1x1-1))] 

= 2tr [d>|x|-i] -tr [d>(|x|-i-Z^,(|x|-i))] . 

According to [50j (see also (|102p ') the first term on the right side is 0{Z'^/^). More- 
over the second term is ©(Z^i/i^) by Lemma hence much smaller than the first 
term. Now, 

(61) D{py - pu^>,P> + Pu.>) = ^tr[dy {v~Uc{v))] , 

Decomposing the trace in (I6ip into the orbitals contributing to (i> and scaling 
X — !■ x/c enables us to employ Lemma 12.91 to obtain the bound 

tr [d> (v - Uciv))] < ^^«(0) tr [d>p2] . 

This concludes the proof, since from |5D] (see (|102p ) we conclude that the trace on 
the right-hand side is 0(Z'^/3). □ 
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4.2.2. Contribution from low angular momenta to the Coulomb energy. We now 
show that the term TZs in (f60|) is neghgible. 

Lemma 4.8. In the limit Z — s- oo one has uniformly in k ^ Z/c G (0, k^] 

Proof. We first treat the term D{p^ ^, pjj y + p>). By construction the densities 
p^ j are sphericaUy symmetric and satisfy according to (|59p 



(62) J pfji{ii)dx = J pf{x)dx = 2{2l + l){K-l), 0<1<L. 
RecalUng the choice of K and L we see that 

(63) y'p^,<(x)dx = 0(Zi/2). 

It foUows from p02p and Lemma [TBI that 

J |x| 
Hence Newton's theorem [44] yields 

D{Pu,<,Pu,> + P>) < ^ J Pc/,<wdx y j^;^ dy = C'(Z 

In the remainder of the proof we are concerned with the term D{p^ ^, p^ ^ + P<)- 
Noting that 

D{p§,<,P§,<+p'<) < ^^(P^,<,P^,<) + ^^(P<,P<)- 

and that according to [SUl Prop. 3.5] D{p^,p^) = 0{Z^^^^), it suffices to consider 
D{Pu <' Pu <)■ We spht the lowest angular momentum corresponding to ^ < 2Z/c— 
1/4 =: ?o off and define 

and 

:= ti-c^ (Zic(d?)(x,x)) , pg^H := trc2 (Z^,(d^)(x, x)) . 

Note that in case Iq < there is no need for this procedure. Accordingly, we 
estimate 

^(P^,<,P^.<) <2i?(p^,H,P^,h) + 2i?(pg,H,p^,H). 

For an estimate of the second part corresponding to Iq < I < L, we apply the 
following angular momentum barrier inequality 

(64) i3,(0)>Z^e(^^X{|x|<n}) 

on where r/ = {{I + Xj^fc^ - AZ'^) /{AZc^) and I > 2Z/c. This bound follows 
by applying Z^i to the inequality in [27l Lemma 2.6] with Ri = [{1 + 1/2)'^-4k'^]/{Ak) 
and scaling x i— > x/c. 
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Inequality ([64]) implies 



1 



tr [U, (|xr 1) df] < — tr [BM df] + tr [U, (|x|-ix{|x|>n}) di 
1 iZ 

Here the last inequality used the fact that eigenfunctions of df are eigenf unctions 
of i3c[^|x|~^] with negative eigenvalue. Now, note that 

{I + 1/2)2 _ iz'^/c^ = {1 + 1/2 + 2Z/c){l + 1/2 - 2Z/c) > const {I + l/2f 

since I > Iq. Hence, using (|62[) and summing over I we obtain 

f ^^ij^dx = |] tr [Z^c(|xri)df] < const + 1/2)-' / pf (x)dx 

Accordingly, Newton's theorem and ([55]) yield 

^(P^H^P^h) < ^ / P^,H(x)dx I ^^dx = 0(Z^V6iogZ). 

Finally, we consider the contribution from I < Iq. Note that then I < 2k^ — 1/4 < 
2. We claim that the electrostatic energy corresponding to the electrons in this 
subspace is bounded by 

(65) D{p§^^,pfj^^) < const cK^. 

Since by the choice of Iq one has 2Z/c > / + l/4> 1/4, estimate (|55|) will imply that 
D{pf \-,Pij I-) < const ZK^ = ©(Z^/'^) and hence complete the proof of Lemma 
14.81 By scaling it suffices to prove (|65l) for c = 1, which we will assume in the 
sequel. The Hardy-Littlewood-Sobolev inequality (see, e.g., [351 Thm. 4.3]) implies 
that 

(66) < const ||p^,h||6/5- 

The triangle inequality together with the definition of U and ((74|) yields 

(67) l|PC/.h||6/5< J2 E ll*-^"ll?2/5' 

where {■(/'q|Q! G A} stands for the collection of normalized eigenfunctions building 
up df , i.e., the corresponding sum ranges over all indices {j,l,m,n). We further 
estimate with the help of Lemma IB. II and Theorem 12.41 

\\'^u^a\\l2/5 ^ const IIV'allL/S ^ COnst . 

This, together with ([55]) . (|67p and the fact that the number of indices in A is 
bounded by a constant times K proves ([65]) . □ 



SCOTT CORRECTION 



29 



4.2.3. Finishing the proof. We repeat ([60]) 

E^iZ) - E^iZ) > tr[5[^]4] ~ tr[Be[||]d^] - 7^l - 7^2 - 7^3 - 7^4 - const 
By Lemmata 14.61 H771 and 14.81 we have uniformly in k ^ Z/c E (0, k^] 

7^l = 0(^23/12)^ 7^2 = 0(^5/^), 7^3 = o(z"/^iogZ), 

so these terms are of lower order than Z^^/24^ Next we scale x i-^ x/c and obtain 

tr[5[Z|x|-i]4] - tr[Bc[Z|x|-i](if ] = Z2s(k) - 7^4 
where s(k) is introduced in ([8]) and 

L — 1 oo 

7^4 :-Z2«;-2^(2? + l) ^ ^ (A„(s,(ac))-A„(6,,(«))- 

1=0 j=l±l/2n=K-l 

By Theorem 12.11 there is a constant such that for all < k < 

L — 1 oo 

o<7^4<z2K-2^(2^ + l) ^ ^ |A„(6,v('«)l 

/^O j^l±l/2 n^K~l 

L— 1 oo 

<constZ2^(2/ + l) {n + iy^ 

l=Q n=K-l 

< const Z^L^K-^ = ©(Z"/^). 
This concludes the proof of the lower bound and hence of our main result. □ 
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Appendix A. Partial wave analysis 

For the convenience of the reader and for normalization of the notation we gather 
some fact on the partial wave analysis of the Brown-Ravenhall operator. 

We denote by Yi^m the normalized spherical harmonics on the unit sphere §2 
(see, e.g., [43], p. 421) with the convention that Yi^m = if |m| > I, and we define 
for j G No + 5, Z G No, and m — ~j, . . . , j the spherical spinors 

/ /.7+m V I. .\\ 
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The set of admissible indices is X :— {{j,l,m) \j e N — 1/2, I = j ± 1/2, m = 
—J, j}. It is known that the functions fljd,m, (j, I, tti) G T, form an orthonormal 
basis of the Hilbert space L^(S^;C^). They are joint eigenfunctions of J^, J3, 
and with eigenvalues given by j(j + 1), Z(/ + 1), and m. The subspace ^j,i,m 
corresponding to the joint eigenspace of total angular momentum with eigenvalue 
j{j + 1) and angular momentum with eigenvalue 1(1 + 1) is then given by 

= span{x ^ |x|"i /(|x|) %,/,,„(wx) | / e 

where Wx := x/|x|. This leads to the orthogonal decomposition 

j 

(69) -^=0 %^ %;=0%i,m, 

of the Hilbert space of two spinors. 
We note that the Fourier transform, 

(70) t/'(p) (27r)~3/2 /■ e-'P'^^Cx) dx, 

leaves the spaces Sjj^i invariant. Namely, if we decompose ip according to (pS]) . 

then 

with the Fourier-Bessel transform 

(72) ^r'^l^ /(r)j,(rp)rpdr. 
Here j/ is a spherical Bessel function. Moreover, 

poo nOC 

(j,/,m)el ° (i,!,m)el ° 

Appendix B. Properties of the twisting operators 
We define the helicity operator H = ojp • cr on ^ by 

(73) HVXp) :=<T.u;p^(p). 
It follows from the pointwise identity 

(74) (Wp • Cr)fljj.m{uJp) = -rij,2j-i,m(Wp), 

see, e.g., Greiner [29l p. 171, (12)], that H is an isomorphism between and 
f)j,2j-i- Moreover, since [cr ■ a) [cr ■ b) = a ■ b + icr ■ [a y. b) for any a,b ^ M.'^, we 
infer that H is an involution on Sj, i.e., H = H^^. 

We shall need to consider H on spaces with p ^ 2. The relevant properties 
are summarized in the next lemma, together with those of the operators 

(75) C^(p) :=<i>,(p)V'(p), 

introduced in Note that while $0 acts trivially on the spin, $1 involves the 
helicity H. 
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Lemma B.l (i^'-properties of H and ^^). The operators H and ly ~Q,1, 
extend to bounded operators from LP{M!^ , C^) to LP{M!^ , C^) for any p € (1, oo). 

Proof. The LP-boundedness of H follows from that of the Riesz transformation, 
see [6TJ Ch. II-III]. Therefore, to prove the statement about the operators it 
suffices to consider the operators (j>i, defined analogously as in ([75| on L^(M'^). Since 
p I— > (j>iy{p) is smooth away from the origin and p^d^(f>i, is bounded for fc = 0, 1, 2, 
the Hormander-Mihlin multiplier theorem [611 Thm. IV. 3] implies that (j)i, extend 
to bounded operators from LP(R'^) to Lp(R^) for any p e (1, oo). □ 

Lemma B.2. For all p, q e 

(76) l-$o(p)*o(q) -*i(p)*i(q) 

= \jz ('J'-(p) - 'J'-(q))' + \ (*i(q)*i(p) - *i(p)*i(q)) • 



and furthermore 



v=0 



|$o(p)-*o(q)l'< "^'^ 



2 



l^l(q)<5>l(p)-a>,(p)^,(q)|<^|||-^ 

Proof. The first equality is an immediate consequence of the definition of $o and 
$1. From this definition we also conclude by an explicit calculation that 

(77) |$o(p)-$o(q)r = (0o(p)-<^o(q))' < 'P""^'' 



8i?(p)2£;(q)2- 
Moreover, for a proof of the next inequality we write 

|$i(p) - $i(q)|2 ^ (0,(p) _ 0,(q))2 + 0i(p)0i(q)|,^p _ 

and estimate the last two terms with the help of the inequalities 

(78) (0i(p)-0i(q)) < 2E{prEiqr " 2£;(p)2i?(q)2 ' 

and 

Finally, for a proof of the last inequality we use 

|$i(q)$i(p) - $i(p)$i(q)| = 2(/.i(p)0i(q) |cr • {ujp X tjq)| < 2,^i (p)(/.i (q)|wp-Wq| . 
Using again ([79]) concludes the proof of the third inequality. □ 
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Appendix C. Basics of relativistic hydrogenic operators 



In this section we collect - following [M] - some basic properties of the operators 
-Bfj and which describe hydrogenic atoms in the Brown-Ravenhall respectively 
Chandrasekhar model. For pedagogical reasons we first discuss their massless ana- 
logues, 



(80) IpI - 5 (Ixl-i 



C.l. Massless case. Expanding ip as in ([7T|) and using l|74p yields [2] the follow- 
ing partial diagonalization of the massless operators, 

Here the operators &^°''(k) and cf'\K) are densely defined in L^(R+) through their 
quadratic forms. 



(/,c[°H«)/> 



p|/(p)Pdp- K 
p|/(p)Pdp- K 



OO /'OO 



JO 

OO poo 



f{q)kj {q,p) f{p)dqdp, 



f{<l)h {q,p) fip)dqdp, 



with maximal form domain denoted by £l{bj^\K)) and !Cl{cf \K)). In the above 
expression, the integral kernels and kf are given by 



(81) 



kfip,q) 



+ 1/2 I 2 



1 / P 



(82) fcF(p,'?):-^Q^(^(f + f)), 

where Qi are Legendre functions of the second kind, i.e.. 



(83) 



Pi{t){z~t)-^dt 



with Pi standing for Legendre polynomials; see Stegun [5D] for the notation and 
some properties of these special functions. 

It was proven in [14] and [35l Eq. (5.33)] that ([80)1 are self-adjoint and lower 
bounded if and only if k < k^, # = -S,C, cf. p^ . More can be said about 
the reduced operators 6^° "*(«;) and c[°^(k). They are lower bounded (in fact, non- 
negative) if and only if 

(84) 
(85) 



t^j Jo 



1 



dt 



This follows by the same lines of reasoning as in [14]. 
Since [60l (8.4)] Pn{t) ^l,P^{t)^ t, we have 

(86) Qo(t) = ^log^, Qi(t) = ilog^-l, 



t-1 
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such that K^' ^ 2/7r, = 7r/2 and thus ^ 2/ {2/tt + tt/2). 

The critical coupling constants and are strictly increasing in j and / and, in 
particular, ~ ^^"^ ~ ■ This follows from the pointwise monotonicity 

(87) Qi{t)>Qi'{t) for > ? and t > 1 

which, in turn, is evident from the integral representation 



Qi{x) ^ j — ===^=dz, x>l] 

Jx+V^^^ V 1 - 2xz + 

see Whittaker and Watson [TOl p. 334, Chap. X, Sec. 3.2]. 

C.2. Massive case. Similarly as before one obtains the following partial diagonal- 
ization of the massive hydrogenic Brown- Ravenhall and Chandrasekhar operators, 

(88) (V', i?KV')= X! {^l'^j,m,l,bj,l{K) J^l^j^rri.l), 

(L7n,s)el 

(89) {iJj^Ci^iIj) = ^ {J^l1pj,m,l,Cl{K)Tlll)j,m,l)- 

Here the operators bj^i{K) and q(k) are densely defined in L^(M_|_) through their 
quadratic forms, 



(90) if, bMf) -.^ {E{p)-l)\f{p)\'dp-K / f{q)k"{q.P)f{p)MAp. 



(91) {f,ci{K)f):= {E{p)~l)\f{p)\'dp-K / f{q)k't{q,p)f{p)dqdp 

JQ JQ Jo 

with maximal form domain denoted by £l{bj^i{K)) and 0(q(k)), cf. |14| . In the 
above expression, the integral kernel k^i depends, in contrast to the massless case, 
on both j, I, and is given by 

kf,{p,q) i [mp)Qi {k + f)) Mq) + Mp)Q2,-i + f)) ^1(9)^ ■ 

The form (|90|) defines a self-adjoint semi-bounded operator bj^i{K), if and only if 
K < nf (Evans et al. [H]). In fact bj.i + is positive (Tix [69 ). A trivially 
modified argument shows that ([91]) defines a self-adjoint semi-bounded operator 
ci{k), if and only if k < . 

In fact the semiboundedness of the massive cases and the massless cases are 
equivalent, since the differences of the massive and massless forms are bounded 
(Tix [IHl Thm. 1]). 

Appendix D. Critical Chandrasekhar operator on a finite domain 

Lieb and Yau [41] have shown that the critical Chandrasekhar operator |p| — 
k'^|x|~^ when restricted to a ball has only discrete spectrum with eigenvalues ac- 
cumulating at infinity at the rate predicted by the semiclassical result for |p| alone. 
This is remarkable since the semiclassical phase-space volume corresponding to 
IpI — k|x|~-'^ is infinite. 

We aim to prove an analogous result for the Chandrasekhar operator restricted 
to the fixed angular momentum subspace corresponding to Z = 1 and finite domain. 
In the proof of Theorem 12.11 it is essential to handle coupling constants which are 
larger than n'^ , all the way up to and including k^. 
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In order to define the above operator we consider for i? > and I G N tlie Hilbert 
space 

UR) {/ e L^O, c») I (^fV) (r) = for all r > R} , 
where JF/ denotes the Fourier-Bessel transformation, cf. (|72p . The quadratic form 
given by {f,cf^\K)f) with domain di{R) n0(c|°^(K)) defines for all k < Kp a 
self-adjoint, non-negative operator in which we will denote by cf^\K,R). 

Lemma D.l. Let I G N. There exists some constant such that for all R > 0, 
fi > 0, and K < Kp 

(92) tr (c^°^ (k, R) - /i^ < const R. 

We have not tried to track the ^-dependence of the constant, since the cases 
1 — 0,1 will be enough for our purpose. 



Proof. For a proof of ([921) we basically follow the argument in [41]. The starting 
point is the following reduction to a simpler variational problem involving only 
functions. Namely, for any non-negative function h : K.+ —^ K+, let 

C l-oc 

Then 



(93) 



tr (q(kP, R) - > inf |^ a{p) {p - - t{p)) dp | < cr < Af;| 



where Mi :— R sup^yQ{2 / TT)r'^ jf (r) . The proof of ([93]) is analogous to the one of 
[4H Eq. (7.8)]. We merely replace the Fourier transformation in M.'^ by the Fourier- 
Bessel transformation !Fi in M+ . 

From now on we assume that / > 1 and comment on the necessary changes in 
case / = at the end. We choose h of the form 

p-i-(A/2)p-2 [ip>A, 
(2A)-i iip<A. 

Below we shall show that the constant A can be picked in such a way that for some 
(5 > 

ifp>d-^A, 
-const if p < 6^^ A. 

In view of (|93p this will prove the result, since then 



h{p) = 



(94) p-fi~ tip) > 



inf I J a{p) {p~ ti- t{p)) dp I < < Af; I > -const A^^^i^Mi j dp 



-const 5^^ I? Ml 



To prove ([94]) we recall that Qi{\{t + })) f = 7i"(Kp)-i, cf. ([SS]), and hence 
by a straightforward calculation 



dt 



piT^pW)-F{A/p)) ifp>A, 
2\{-F{l) + Fip/A)) ifp<A. 
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Here for < s < 1 we have set 

F(s) :^ rQ,(i(t+i))(i-i)^dt. 

Since Qi{t) < Qi{t), which vanishes Hke a constant times as r ^ oo, one has 
F{s) ^ as s ^ 0. Choosing S S (0, 1) such that F{s) < ^F{1) for aU < s < 5, 
we have shown that for all p > S^^A one has 

A F(l) , F(l) 

For A < p < S^^A we use the monotonicity, dF/ds > 0, to bound 

p - t{p) > 0. 

FinaUy, for < p < A we drop the term F{p/A) > to obtain 

P't{p)>~^F{l)>-A^^^^ 



2 ' ' - 2 
Choosing A :— 8/i/i^(l) yields the claimed inequality ((94|) . 

In case / = 0, the function h can be chosen as before. However, the corresponding 
expressions F{1) — F{s) should be interpreted as a single integral, and estimated 
with slightly more care. □ 

Corollary D.2. Let I G N. Then there exists some constant such that for all 
< K < nf, all n > and all functions x on R-|_ which satisfy % > on [0, R) and 
X = on [R, oo) for some R > 0, one has: 

NiiO, X (IpI - Klxj-i - /i) x) < const fiR. 

Proof. The variational principle implies that 

Ni{Q, X (IpI - «:|xr' -i^)x)< Ni^i, c[°)(«, r)). 

Indeed, if V; is the negative spectral subspace of x (|p| ^ '^1^1^^ ^ y") X with fixed 

I, then any / e J'lxVi C Si{R) satisfies (/, {c\°\k, R) - /i) /) < 0. 
The assertion now follows from 

iV(^,c[°^(K,i?)) < const fiR. 

For a proof, we note that the elementary inequality X(-oo,/j)(^) ^ ^^x^^^ ' valid 
for any /i < A, together with Lemma FD . 1 1 implies that 

N{fi,c\°\K,R)) < {X - fj.)-HT{cf\K, R) - A)_ < const (A - /z)"^ A^i?. 
The proof is completed by optimizing over A. □ 

Appendix E. The trial density matrix 

In this section we define the density matrices d^ and d^ that we use to bound the 
Schrodinger energy, respectively the Brown-Ravenhall energy, from above. Both 
density matrices are split into two parts corresponding to low and high angular 
momenta 

d''^ := d%+d>, d^ -.^d^+dy. 

Low angular momenta correspond to orbits whose perinucleon is close to the nu- 
cleus, while high angular momenta ensure that the orbits are never close to the 
nucleus. We will cut between these two at L := [Z^^^^]. 
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E.l. Low angular momenta. In the vicinity of the nucleus the nuclear attraction 
dominates the interaction with the other electrons. This motivates to choose the 
orbitals as the ones of the Bohr atom, i.e., as the eigenfunctions of the unscreened 
operator with nuclear charge Z . The corresponding density matrices c?< are of the 
form 

^-E^f' 4= E 4 

i=0 J=i±l/2, J>l/2 



and 



E E 



m— — j n— 1 

Here K = [const Z^^"^] with some positive constant, i.e., on the order of the last 
occupied shell of the Bohr atom. We now turn to the definition of the orbitals 
^ „ for which we consider the cases ^ = B,S separately. 
In the Brown- Ravenhall case we choose tp^i ^ „ such that its Fourier transform 

is 

where ffin is the n-the eigenfunction of the operator Vcbjj{Z/c)V* in L^(R_(_). 
Here the unitary scaling operator Vc is defined by {Vcf){p) :~ cr^/"^ f{p/c) and 
we recall that the operator bjj{K) was defined in Subsection IC.2I The operators 
Vc bj,i{Z/c) V* appear as the angular momentum reductions of i?c[^|x|~^]. Indeed, 
by ([88| and scaling one has 



B,[Z|x|-l] V) = E (V-J-.™.'' b,,l{Z/c) V: i^j,,n,l). 



In the Schrodinger case we choose 

where ff^ is the n-th eigenfunction of —^372 + ^^^77^ ^ "f in L^(^+) with Dirichlet 
boundary conditions. 

E.2. High angular momenta. For large angular momenta, the electrons are suf- 
ficiently far from the center moving ~ classically speaking - slowly. This motivates 
to pick non-relativistic orbitals in both in the relativistic and non-relativistic case. 
Moreover, for large quantum numbers the correspondence principle would predict 
quasi-classical behavior (in the quantum sense) as well. This motivates the following 
choice which we take - with slight modifications - from [50] : 

j 

(95) (i>:=E'^'' •= E E Wn,l\fn,l^j,l,m){Vn,l^j,l,m\- 

1>L j=l±l/2rn=~j n£Z 

We repeat at this point the construction of the Macke orbitals (pn^i and their weights 
Wn,i- We will also present a new estimate not directly given in that paper. 

The semi-classical mean-field in which the electrons move is the Thomas-Fermi 
potential (/>tf (see (|49)) V According to Hellmann ^32] the semi-classical electron 
density for fixed angular momentum is 



(96) af(r) 




nz(pTF[r) r-^ — 



SCOTT CORRECTION 



37 



where we added the factor nz ~ {1 — aZ^^^'^^l'^ for normahzation purposes with 
some fixed positive a and where we replaced the self-generated field of the sum of 
the radial densities cr; by the Thomas-Fermi potential. We will write for the 
functions when a = 0, i.e., no normalization factor occurs. In passing we note 
that the densities are the minimizers of the Hellmann functional with external 
potential given by the Thomas-Fermi density and no other interaction between the 
electrons (see [56 ). 

The functions vanish for large I and we define 

h! min{? G N | cif = 0}. 

By scaling, fc' is of the order Z^l'^ . Moreover, since the function r ^ (f>TF{r)r^ has 
exactly one maximum, the support of is an interval [ri(/), r2(/)]. 

We cannot use the density aj^ directly in defining semi-classical orbitals, since 
the derivative of its square root is not square integrable. Thus we pick two points, 

(97) xi{l) := nil) + T{1 + i)Z-\ X2il) := r^il) - SZ'^'^ 
for some positive S and T G (0,4), and set 

r2(2; + l)aV+2, re[0,xi(0], 

(98) P;(r) := L/^(r), r e [x^{l) , x^il)], 

[2(2/-Hl)/32exp(-23/2z2/3r), re [x2{l),oo). 

The constants a and (3 are chosen such that pi is continuous. We suppress their 
dependence on I in the notation. 

Next, we define for I < k' and ri e Z the Macke orbitals 



'^-'-^=2(2rM)./o 



(99) (p„,(r) := Va«^i.fc„,c,W 

r 

where Q '■ [0, oo) — > [0, 1) is the Macke transform 

(100) Ci{r) oo 

Jo Pi(i)di 

For l> k' we set ipn,i ■= 0. The integral 

/>oo 

2(2r 

which is independent of j and m, will represent the number of electrons in the 
angular momentum channel {j,l,m). Moreover, we set ei := Nj^i^m — [Nj^i^m]- If 
[Nj.i.rn] is odd, we pick kn^i = 2n, otherwise = 2n — 1. The weights are chosen 
as 

fl \kn.l\<[Nj,Lm]~l 

(101) Wn,i := I Si/2 \kn,i\ = [Nj,i^m] + 1 

[ otherwise 

which guarantees that X)nez ~ ^j-J,m- 

Strictly speaking, our trial density matrix differs from the one used in [50 , since 
we label the orbitals by the modulus of total angular momentum, by the third 
component of total angular momentum, and by the orbital angular momentum. 
This, however, is merely a minor rearrangement of terms. 
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We also adapt to atomic units used in this paper which changes the value of the 
Thomas- Fermi constant and gives a factor 1/2 in front of all three kinetic energy 
terms in the Hellmann-Weizsacker functional. 



E.3. Energy estimates for high angular momenta. For the convenience of the 
reader, we gather from [5(T (based on the construction in 55 ) two estimates on 
the order of the average kinetic and potential energy of the Schrodinger operator 
associated with the semi-classical density matrix (i> , 

(102) tr(p2d>) =0(Z^/3), tv{\^\-^dy) ^0{Z^/^). 

We also need a more detailed estimate on the kinetic energy. 

Lemma E.l. Let L = [Z^/^'^]. Then for large Z, 



(103) 



J2l'^^<P^di) = OiZyL). 



l=L 



Proof. The definition of di implies (cf. [50l (2.3)]) that for angular momenta I < k' 
one has 



(104) tvip^di) 
where we set 



,2 a; 3 , ^(^ + 1) 



dr + Fi 



„ oil 
Fi := — 
3 



iV2 

j,l,m 



-1 + 6ei - 3ef _^ 2ef - 6ef + 4e/ 



Pfdr, 



and emphasize that ai should not be confused with a from 
Proposition 3.6] 



ai : 



4(2/-!- 1)2' 
According to [SHI 



l^^ Fi<^Fi< const Z 



5/3 



l=L 



l=L 



where L — [Z^^^^]. The first term on the right-hand side of (|104p is estimated 
according to 



(105) 
with 



ai 3 /(/ + 1) 



-Pi 



-Pi 



dr 



< 



dr + Gi+Hi+Ii. 



Gi := 



H, 



X2(l) 



'2 I 

+ ^Pl 



-pi 



dr < const Z^ I 



-3/2 



dr < const Z"^/^ ' ' + 2 
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where the inequahties were obtained by integration as in [501 (3.4)]. Inequahty [501 
(3.9)] for the gradient term in the middle region reads 



h / ^^Pl ' dr < const M + - 



-1/2 



2/ V 2/ min{Z + i,Zi/4} 



This imphes 

^ r"^ Gi < const < Z^L-^/^, 

oo k 

l^^ Hi<Y Hi < const Z'^/s log k' < const Z'^/^ log Z 

l=L l=L 
oo 

Y h < const + Z log Z + ^2^-5/2 + z5/3L-3/2j < const Z^^l^\ 

l=L 

It thus remains to estimate the sum of the first terms on the right-hand side of 
(|105p . We begin with the first summand, 

-1 /-oo -1 />oo 

E4/ fpf(r)Mr< const (Z/r - Z2/,2)3/2^, 

1 /-oo 

= const ^"'^'(1 - r-i)'/'dr = ©(Z^/L) 



where we used that the Thomas-Fermi potential is bounded from above by Z /r. 
This leaves the second summand, 

go -, nOQ (1 I ]_\2 oo „oo 

^ ^ / ^ pf(r)dr < const ^ ? / r-2(Z/r - l^/r^ fl^dr 



-.L l=L 

-5/2n _r-i^i/2 1^ ^ .orz2, 



-constZ^VA / r-^/2(i_r-i)V-dr = 0(Z2/L), 

which completes the proof of Lemma IE. II □ 

References 

[1] Volkcr Bach. Error bound for the Hartree-Fock energy of atoms and molecules. Comm. Math. 

Phys., 147:527-548, 1992. 
[2] A. A. Balinsky and W. D. Evans. On the virial theorem for the relativistic operator of Brown 

and Ravenhall, and the absence of embedded eigenvalues. Lett. Math. Phys., 44(3):233— 248, 

1998. 

[3] Hans A. Bethe and Edwin E. Salpeter. Quantum mechanics of one- and two-electron atoms. 
In S. Fliigge, editor, Handbuch der Physik, XXXV, pages 88—436. Springer, Berlin, 1 edition, 
1957. 

[4] G. E. Brown and D. G. Ravenhall. On the interaction of two electrons. Proc. Roy. Soc. 
London Ser. A., 208:552-559, 1951. 

[5] Roch Cassanas and Heinz Siedentop. The ground-state energy of heavy atoms accord- 
ing to Brown and Ravenhall: absence of relativistic effects in leading order. J. Phys. A, 
39(33):10405-10414, 2006. 

[6] Charles G. Darwin. The wave equation of the electron. Proc. Roy. Soc. (London), 1928. 



40 R. FRANK, H. SIEDENTOP, AND S. WARZEL 

[7] Ingrid Daubcchics. An uncertainty principle for Fermions with generalized kinetic energy. 

Comm. Math. Phys., 90:511-520, September 1983. 
[8] .lean Dolbeault, Maria J. Esteban, and Eric Sere. On the eigenvalues of operators with gaps. 

Apphcation to Dirac operators. J. Fund. Anal, 174(l):208-226, 2000. 
[9] Berthold-Georg Englert. Semiclassical Theory of Atoms. Number 300 in Lecture Notes in 
Physics. Springer-Verlag, Berlin, 1 edition, 1988. 
[10] Berthold-Georg Englert and Julian Schwinger. Statistical atom: Handling the strongly bound 

electrons. Phys. Rev. A, 29(5):2331-2338, 1984. 
[11] Berthold-Georg Englert and Julian Schwinger. Statistical atom: Some quantum improve- 
ments. Phys. Rev. A, 29(5):2339-2352, 1984. 
[12] Berthold-Georg Englert and Julian Schwinger. Atomic-binding-energy oscillations. Phys. Rev. 
A, 32:47-63, 1985. 

[13] A. Erdelyi, W. Magnus, F. Oberhettinger, and F. G. Tricomi. Tables of integral transforms. 
Vol. II. McGraw-Hill Book Company, Inc., New York- Toronto-London, 1954. Based, in part, 
on notes left by Harry Bateman. 

[14] William Desmond Evans, Peter Perry, and Heinz Siedentop. The spectrum of relativistic one- 
electron atoms according to Bethe and Salpeter. Comm. Math. Phys., 178(3):733-746, July 
1996. 

[15] C. FefFerman and L. Seco. Eigenfunctions and eigenvalues of ordinary differential operators. 

Adv. Math., 95{2):145-305, October 1992. 
[16] C. Fefferman and L. Seco. The density of a one-diincnsional potential. Adv. Math., 107(2):187— 
364, September 1994. 

[17] C. Fefferman and L. Seco. The eigenvalue sum of a one-dimensional potential. Adv. Math., 

108(2):263-335, October 1994. 
[18] C. Fefferman and L. Seco. On the Dirac and Schwinger corrections to the ground-state energy 

of an atom. Adv. Math., 107(1):1-188, August 1994. 
[19] C. Fefferman and L. Seco. The density in a three-dimensional radial potential. Adv. Math., 

111(1):88-161, March 1995. 
[20] C. L. Fefferman and L. A. Seco. An upper bound for the number of electrons in a large ion. 

Proc. Nat. Acad. Sci. USA, 86:3464-3465, 1989. 
[21] C. L. Fefferman and L. A. Seco. Asymptotic neutrality of large ions. Comm. Math. Phys., 

128:109-130, 1990. 

[22] C. L. Fefferman and L. A. Seco. On the energy of a large atom. Bull. AMS, 23(2):525-530, 
October 1990. 

[23] Charles L. Fefferman and Luis A. Seco. Aperiodicity of the Hamiltonian flow in the Thomas- 
Fermi potential. Revista Mathemdtica Iberoamericana, 9(3):409-551, 1993. 

[24] E. Fermi. Un metodo statistico per la determinazione di alcune proprieta dell'atomo. Atti della 
Reale Accademia Nazionale dei Lincei, Rendiconti, Classe di Scienze Fisiche, Matematiche 

e Naturali, 6(12):602-607, 1927. 
[25] E. Fermi. Eiiic statistische Begriindung zur Bestimmung ciniger Eigenschaften des Atoms 
und ihre Anwendungen auf die Theorie des periodischen Systems der Elemente. Z. Phys., 
48:73-79, 1928. 

[26] Rupert L. Frank, Elliott H. Lieb, and Robert Seiringer. Hardy-Lieb-Thirring inequalities for 
fractional Schrodinger operators. Joum. American Math. Sac, Online, 2008. 

[27] Rupert L. Frank, Heinz Siedentop, and Simone Warzel. The ground state energy of heavy 
atoms: relativistic lowering of the leading energy correction. Comm. Math. Phys., 278(2):549— 
566, 2008. 

[28] Walter Gordon. Die Encrgiciiiveaus des Wasserstoffatoms nach der Diracschen Qaunteiitlic- 

orie. Z. Phys., 48:11-14, 1928. 
[29] Walter Greiner. Relativistic Quantum Mechanics, volume 3 of Theoretical Physics - Text and 

Excercise Books. Springer, Berlin, 1 edition, 1990. 
[30] Marcel Griesemer, Roger T. Lewis, and Heinz Siedentop. A minimax principle for eigenvalues 

in spectral gaps: Dirac operators with Coulomb potential. Doc. Math., 4:275-283, 1999. 
[31] Marcel Griesemer and Heinz Siedentop. A minimax principle for the eigenvalues in spectral 

gaps. J. London Math. Soc. (2), 60(2):490-500, 1999. 
[32] Heinrich Hellmann. Ein kombiniertes Naherungsverfahren zur Energieberechnung im Vielelek- 

tronenproblem. H. Acta Physicochim. U.S.S.R., 4:225-244, 1936. 



SCOTT CORRECTION 



41 



[33] Webster Hughes. An Atomic Energy Lower Bound that Gives Scott's Correction. PhD thesis, 

Princeton, Department of Mathematics, 1986. 
[34] Webster Hughes. An atomic lower bound that agrees with Scott's correction. Adv. in Math., 

79:213-270, 1990. 

[35] Tosio Kato. Perturbation Theory for Linear Operators, volume 132 of Grundlehren der math- 
ematischen Wissenschaften. Springer- Verlag, Berlin, 1 edition, 1966. 

[36] W. Lenz. Uber die Awendbarkeit der statistischen Methode auf lonengitter. Z. Phys., 77:713- 
721, 1932. 

[37] Elliott H. Lieb. Variational principle for many-fermion systems. Phys. Rev. Lett., 46(7):457— 
459, 1981. 

[38] Elliott H. Lieb and Michael Loss. Analysis, volume 14 of Graduate Studies in Mathematics. 
American Mathematical Society, Providence, RI, second edition, 2001. 

[39] Elliott H. Lieb and Barry Simon. The Thomas-Fermi theory of atoms, molecules and solids. 
Adv. Math., 23:22-116, 1977. 

[40] Elliott H. Lieb and Walter E. Thirring. Inequalities for the moments of the eigenvalues of 
the Schrodinger Hamiltonian and their relation to Sobolev inequalities. In Elliott H. Lieb, 
Barry Simon, and Arthur S. Wightman, editors. Studies in Mathematical Physics: Essays in 
Honor of Valentine Bargmann. Princeton University Press, Princeton, 1976. 

[41] Elliott H. Lieb and Horng-Tzer Yau. The stability and instability of relativistic matter. 
Comm. Math. Phys., 118:177-213, 1988. 

[42] Paul Mancas, A. M. Klaus Miiller, and Heinz Siedentop. The optimal size of the exchange 
hole and reduction to one-particle Hamiltonians. Theoretical Chemistry Accounts: Theory, 
Computation, and Modeling (Theoretica Chimica Acta), 111(1):49— 53, February 2004. 

[43] Albert Messiah. Mecanique Quantique, volume 1. Dunod, Paris, 2 edition, 1969. 

[44] Isaac Newton. Philosophiae naturalis principia mathematica. Vol. L Harvard University 
Press, Cambridge, Mass., 1972. Reprinting of the third edition (1726) with variant read- 
ings. Assembled and edited by Alexandre Koyre and I. Bernard Cohen with the assistance of 
Anne Whitman. 

[45] Thomas 0stergaard S0rensen. Towards a Relativistic Scott Correction. PhD thesis, Aar- 
hus Universitet, http:/ /www.imf.au.dk/cgi-bin/viewers/viewpublications.cgi?id=79, October 
1998. 

[46] Thomas 0stergaard S0rensen. The large-Z behavior of pseudorelativistic atoms. J. Math. 

Phys., 46(5):052307, 24, 2005. 
[47] Julian Schwinger. Thomas-Fermi model: The leading correction. Phys. Rev. A, 22(5):1827— 

1832, 1980. 

[48] Julian Schwinger. Thomas-Fermi model: The second correction. Phys. Rev. A, 24{5):2353— 
2361, 1981. 

[49] J. M. C. Scott. The binding energy of the Thomas-Fermi atom. Phil. Mag., 43:859-867, 1952. 
[50] Heinz Siedentop and Rudi Weikard. On the leading energy correction for the statistical model 

of the atom: Interacting case. Comm. Math. Phys., 112:471-490, 1987. 
[51] Heinz Siedentop and Rudi Weikard. Upper bound on the ground state energy of atoms that 

proves Scott's conjecture. Phys. Lett. A, 120:341-342, 1987. 
[52] Heinz Siedentop and Rudi Weikard. On the leading energy correction of the statistical atom: 

Lower bound. Europhysics Letters, 6:189-192, 1988. 
[53] Heinz Siedentop and Rudi Weikard. On the leading correction of the Thomas-Fermi model: 

Lower bound - with an appendix by A. M. K. Miiller. Invent. Math., 97:159-193, 1989. 
[54] Heinz Siedentop and Rudi Weikard. A new phase space localization technique with application 

to the sum of negative eigenvalues of Schrodinger operators. Annales Scientifiques de I'Ecole 

Normale Superieure, 24(2):215-225, 1991. 
[55] Heinz K. H. Siedentop. On the relation between the Hellmann energy functional and the 

ground state energy of an N-Fermion system. Z. Phys. A, 302:213—218, 1981. 
[56] Heinz K. H. Siedentop and Rudi Weikard. On the leading energy correction for the statistical 

model of the atom: Non-interacting case. Abh. Braunschweig. Wiss. Ges., 38:145—158, 1986. 
[57] B. Simon. Fifteen problems in mathematical physics. In Perspectives in Mathematics. 

Birkhauser, 1984. 

[58] Barry Simon. Functional Integration and Quantum Physics. Academic Press Inc. [Harcourt 
Brace Jovanovich Publishers], New York, 1979. 



42 



R. FRANK, H. SIEDENTOP, AND S. WARZEL 



[59] Jan Philip Solovcj. The relativistic Scott correction. In Maria J. Estcban, Claude Le Bris, and 
Gustavo Scuseria, editors, Mathematical and Numerical Aspects of Quantum Chemistry Prob- 
lems, volume 47/2006 of Oberwolfach Report, pages 52-53. Mathematisches Forschungsinsti- 
tut Oberwolfach, European Mathematical Society, September 2006. 

[60] Irene A. Stegun. Legendre functions. In Milton Abramowitz and Irene A. Stegun, editors, 
Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, 
chapter 8, pages 331-353. Dover Publications, New York, 1965. 

[61] Elias M. Stein. Singular integrals and differentiability properties of functions. Princeton 
Mathematical Series, No. 30. Princeton University Press, Princeton, N.J., 1970. 

[62] J. Sucher. Foundations of the relativistic theory of many-electron atoms. Phys. Rev. A, 
22(2):348-362, August 1980. 

[63] J. Sucher. Foundations of the relativistic theory of many-electron bound states. International 
Journal of Quantum Chemistry, 25:3-21, 1984. 

[64] J. Sucher. Relativistic many-electron Hamiltonians. Phys. Scripta, 36:271—281, 1987. 

[65] Walter Thirring. Lehrbuch der Mathematischen Physik 3: Quantenmechanik von Atomen 
und Molekiilen. Springer- Verlag, Wien, New York, 1 edition, 1979. 

[66] L. H. Thomas. The calculation of atomic fields. Proc. Camb. Phil. Soc, 23:542-548, 1927. 

[67] C. Tix. Lower bound for the ground state energy of the no-pair Hamiltonian. Phys. Lett. B, 
405(3-4):293-296, 1997. 

[68] C. Tix. Self-adjointness and spectral properties of a pseudo-relativistic Hamiltonian due to 

Brown and Ravcnhall. Preprint, mp-arc: 97-441, 1997. 
[69] C. Tix. Strict positivity of a relativistic Hamiltonian due to Brown and Ravenhall. Bull. 

London Math. Soc, 30(3):283-290, 1998. 
[70] E. T. Whittakcr and G. N. Watson. A Course of Modern Analysis; An Introduction to 

the General Theory of Infinite Processes and of Analytic Functions, with an Account of 

the Principal Transcendental Functions. Cambridge University Press, Cambridge, 4 edition, 

1927. 

Department of Mathematics, Princeton University, Princeton, NJ 08544-1000, USA 
E-mail address: rlfrankSmath. prlnceton.edu 

Mathematisches Institut, Ludwig-Maximilians-Universitat Munchen, Theresienstrasse 

39, 80333 Mi;-N(^iiEN, Germany 
E-mail address: h.s@lmu.de 

Department of Mathematics, Princeton University, Princeton, NJ 08544-1000, USA 
E-mail address: swarzelSprinceton.edu 



